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20 .  (continued) 

The  theoretical  formulation  pertains  to  a  bubble  in  a  fluid  which 
Is  Infinite  In  extent  and  contains  no  barriers  or  Interfaces  apart  from 
the  bubble  surface.  An  external  flow  can  be  Imposed,  Imparting  a  stress 
history.  In  addition,  a  spherical  bubble  profile,  obtained  from  present 
experiments  or  theories,  dictates  a  volume  change  and  drives  the  non- 
sphericities  . 

Surface  tension  and  inertial  forces  are  very  important  in  these 
bubbles  and  are  Incorporated  In  the  model.  Because  the  effects  of  fluid 
rheology  are  of  primary  Interest,  the  changes  In  bubble  behavior  which 
occur  between  a  Newtonian  and  non-Newtonian  environment  are  discussed  In 
detail.  An  integral  constitutive  expression  generates  the  stress  field 
in  the  fluid. 

One  concludes  from  the  theory  that  viscosity  and  elasticity  have 
negligible  Influence  on  growth  and  collapse  of  spherical  bubbles  under 
conditions  likely  to  be  encountered  in  practice.  However,  the  behavior  of 
nonspherical  bubbles  in  an  otherwise  quiescent  fluid  is  significantly 
altered  by  the  addition  of  an  elastic  component  to  the  constitutive  rela¬ 
tion,  and  such  changes  are  even  larger  when  an  external  flow  Is  imposed. 

Reproducible  cavities  were  generated  by  optical  cavitation  techniques 
in  both  water  and  dilute  aqueous  polymer  solutions.  Spherical  or  non¬ 
spherical  bubbles  were  produced,  either  near  or  far  from  a  solid  wall. 

Experimental  results  show  good  agreement  with  the  theory  for  non¬ 
spherical  bubbles  In  an  Infinite  quiescent  fluid.  Other  trials,  near  a 
solid  wall,  suggest  the  existence  of  a  competition  between  nonspherical 
modes.  The  Jet  which  is  induced  by  a  solid  wall  near  an  initially  spheri¬ 
cal  bubble  did  not  appear  when  the  bubble  1b  initiated  asymmetrically. 

Finally,  we  have  shown  that  from  a  combination  of  external  flow 
theory  and  experimental  observations  one  can  predict  that  changes  in  rheo¬ 
logical  properties  can  macroscopically  alter  cavitation  behavior  by  chang¬ 
ing  the  influence  of  the  flow  on  the  bubble.  This  is  the  first  work  to 
put  forth  a  plausible  explanation  for  such  an  interaction. 
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I.  INTRODUCTION 

Bubble  dynamics  is  con&only  regarded  as  a  harmless  diversion. 

It  elicits  images  of  soap  bubbles  floating  on  air  currents  and  per¬ 
haps  of  bubble  baths  or  the  gaseous  volumes  which  arise  in  carbonated 
beverages.  However,  bubbles  have  a  myriad  of  beneficial  practical 
applications  as  well  as  potentially  catastrophic  consequences  in  many 
processes  involving  liquid  flow.  Of  particular  relevance  here  is  the 
flow  cavitation  which  can  occur  in  pumping  or  aquatic  propulsion  de¬ 
vices.  Uncontrolled  formation  of  voids  under  these  circumstances  can 
result  in  large  drops  in  efficiency  and  can  eventually  cause  struc¬ 
tural  damage  to  the  apparatus.  Research  activity  has  been  carried  out 
in  this  area  for  most  of  this  century.  Great  strides  have  been  made 
in  reducing  cavitation  damage,  e.g.  the  design  and  uBe  of  super-cavitat- 
ing  propellers  which  induce  void  formation  in  a  controlled  and  predict¬ 
able  manner,  but  understanding  of  these  phenomena  is  far  from  complete. 

A  new  complication  was  introduced  by  the  relatively  recent  avail¬ 
ability  of  synthetic  high  molecular  weight  polymer  molecules  and  the 
subsequent  investigation  of  their  properties.  The  presence  of  "drag 
reducing"  polymers  as  solutes  at  very  low  concentrations  O'-  500  ppm 
levels)  in  aqueous  liquids  dramatically  changes  the  initiation  and 
subsequent  behavior  of  cavitation  in  "water"  tunnel  flows  past  blunt 
bodies.  Attempts  to  extract  the  essential  features  which  determine 
"cavitation  inhibition"  to  allow  some  predictive  capability  for  bubble 
dynamics  and  to  provide  more  fundamental  knowledge  about  the  behavior 
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of  polymer  solutions  in  nonvlscometric  flows  have  achieved  only  limited 
success.  A  new  model  and  experimental  program  is  offered  here  which 
includes  elements  present  in  actual  flow  cavitation*  but  have  not  re¬ 
ceived  much  previous  attention.  For  the  first  time,  dilute  polymer 
solution  non-Newtonian  effects  cause  non-negliglble  changes  in  bubble 
dynamics  which  can,  in  turn,  alter  the  gross  behavior  of  the  void. 

The  features  considered  explicitly  in  this  study  are  presented 
in  Table  1.  Surface  tension  at  the  gas-liquid  interface  is  retained, 
but  heat  and  mass  transfer  effects  are  not  explicitly  included.  Al¬ 
though  Incorporation  of  the  Influence  of  each  of  these  factors  is  not 
original  with  this  work,  the  analysis  of  the  interaction  between  them 
is  unique.  The  simplest  case  of  interest  involves  the  behavior  of  a 
nonspherical  bubble  in  a  fluid  of  infinite  extent  upon  which  no  exter¬ 
nal  flow  is  Imposed.  The  results  show  that  the  oscillations  in  shape, 
which  occur  because  of  the  presence  of  surface  tension,  are  strongly 
coupled  to  the  overall  volume  change  of  the  bubble.  A  change  in  shape 
due  to  fluid  rheology  is  predicted  which  Is  larger  than  any  alteration 
of  the  overall  bubble  size  due  to  the  same  differences  in  fluid  proper¬ 
ties.  The  addition  of  an  externally  imposed  flow  around  the  body,  which 
is  the  only  source  of  asymmetries  in  the  system,  couples  with  fluid 
rheology  to  create  differences  between  model  bubbles  which  are  much 
larger  than  those  changes  without  flow.  The  experimental  investigation 
of  the  behavior  of  an  initially  ncnepherioal  bubble  near  a  solid  wall 
suggests  that  the  well-established  solid  wall  affect  —  the  induction  of 
a  liquid  jet  during  the  collapse  of  an  initially  spherical  bubble  —  is 
reduced  by  the  Initial  asyometry. 
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Features  of  Present  Study 
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These  results.  In  conjunction  with  one  another ,  suggest  s  aech- 
snlsm  by  which  dilute  solution  rheology  Bay  strongly  Influence  flow 
cavitation  behavior.  The  shape  and  subsequent  action  upon  collapse  of 
a  bubble  subject  to  the  velocity  gradients  encountered  during  flow 
cavitation  ere  significantly  changed  by  fluid  rheology.  This  shape  is 
also  important  to  the  interactions  between  the  bubble  and  the  solid 
boundaries  which  are  present.  Thus ,  flow  and  rheology  can  change  gross 
bubble  dynamics.  This  working  hypothesis  is  consistent  with  all  of  the 
present  work  and  also  with  the  puzzling  lack  of  date  showing  rheologi¬ 
cal  effects  in  previous  single  bubble  experiments  and  theories  since 
they  usually  treated  the  factors  Individually. 

In  order  to  demonstrate  the  Importance  of  flow,  asymmetries,  and 
rheology  in  bubble  dynamics,  the  groundwork  laid  by  past  researchers 
Is  very  useful.  From  that  basis,  the  present  aodel  can  be  developed 
with  special  consideration  of  Initiation,  rheological  models  and  the 
methods  employed  to  solve  the  resulting  equations.  Then,  the  means 
by  which  experimental  tests  were  conducted  are  described.  Results  are 
subsequently  presented.  Initially,  they  appear  with  some  consent  as 
to  their  individual  meaning,  then  they  are  discussed  from  a  broader 
perspective.  The  conclusions  so  drawn  are  complete  in  the  context  of 
the  present  work,  but  also  provoke  speculation  and  suggest  future  work. 
Appendices  have  been  included  to  clarify  theoretical,  numerical  and 
experimental  details.  As  a  result,  little  development  of  equations  is 
presented  in  the  primary  text. 
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II.  BACKGROUND 


The  addition  of  small  quantities  of  macromolecules  to  water  (e.g. 
500  ppm  polyacrylamide)  affects  relatively  large  changes  in  the  forma¬ 
tion  and  subsequent  behavior  of  flow-induced  cavitation  bubbles.  This 
phenomenon  has  been  quantified  through  meaaureaent  of  the  cavitation 
inception  parameter,  o which  represents  the  conditions  at  which  cavi¬ 
tation  Initially  occurs  as 


P  -  P, 
s  v 

o,  -  •: - ~5T 


h  P  V2, 


(1) 


where  P  ia  the  static  flud  pressure,  P  is  the  vapor  pressure,  p  is 

8  V 

the  density  and  V  is  the  free  stream  velocity.  Figure  1  shows  typical 
o 

"cavitation  inhibition"  data  for  water  tunnel  flow  past  a  blunt  body. 
The  cavitation  parameter  has  been  reduced  by  as  much  as  70Z  for  poly¬ 
mers  such  as  guar  gum  (Ellis  and  Hoyt,  1968,  Ellis  and  Ting,  1974, 

Oba,  et.  al . ,  1978,  Hoyt,  1976).  Thus,  higher  velocities  or  "more 
severe"  conditions  are  necessary  to  form  voids.  Other  experiments 
show  that  differences  persist  even  after  inception. 

In  water,  the  appearance  of  the  cavitation  bubbles 
is  very  violent  and  chaotic,  consisting  of  many 
very  small  bubbles...  However,  as  the  polymer  is 
added,  the  cavity  looks  more  transparent,  and  shows 
a  regular,  smooth,  wavy  pattern  at  tbs  vapor-liquid 
interface.  (Ting,  1978) 


Gross  differences  appear  in  both  the  nucleatlon  and  subsequent  behav¬ 
ior  of  cavitation  bubbles.  In  this  work,  attention  is  focussed  on  the 
latter  stage,  where  continuum  descriptions  srs  spplled. 


7 


Many  excellent  sources  of  general  information  on  bubbles  ere  avail¬ 
able  Including  volumes  by  Harnmit  and  hia  co-workers  (1980*  Knapp,  st.  at. , 
1970),  as  well  as  review  articles  by  Plesset  and  Prosperetti  (1977), 
Acosta  and  Parkin  (1975),  Plesset  alone  (1977),  and  Prosperetti  alone 
(1981).  All  theoretical  work  is  seen  to  begin  with  Lord  Rayleigh  in  1917. 
He  began  by  assuming  a  spherical  void  to  be  present  In  an  Infinite  medium 
and  neglecting  all  fluid  properties  except  a  constant  density.  A  large 
amount  of  the  bubble  literature  since  that  time  has  been  devoted  to  the 
analysis  of  more  detailed  descriptions  of  fluid  behavior,  including  vis¬ 
cosity  and  surface  tension  and  other  physical  processes,  e.g.  heat  and 
mass  transfer,  compressibility  and  the  transport  of  contaminants.  A 
form  of  these  results  is  employed  in  the  present  work  (Equation  III. 2) 
and  like  most  treatments  neglects  processes  within  the  bubble.  By  re¬ 
taining  an  assumption  of  spherical  symmetry, only  one  spatial  coordinate, 
the  radius  r  in  spherical  coordinates,  has  bearing  on  these  analyses, 
greatly  simplifying  the  equations  of  motion. 

Results  for  spherical  bubbles  are  divided  into  categories.  For 
pure  water,  surface  tension  o  and  viscosity  u  are  weak  functions  of 
other  physical  parameters  such  as  temperature.  This  has  allowed  suc¬ 
cessful  modelling  under  the  assumption  of  constant  values  for  those  two 
coefficients.  However,  the  composition  and  behavior  of  the  pressure  in 
the  interior  of  the  bubble  with  the  evolution  of  the  bubble  is  not  so 
simple.  Plesset  (1977)  has  delineated  two  categories  of  bubble  dynamics 
which  he  has  labelled  gae  bubbles  and  vapor  bubbles.  Gas  bubbles  are 
those  cavities  for  which  the  medium  in  the  interior  is  largely  or  com¬ 
pletely  a  permanent,  noncond ana able  gas.  Por  vapor  bubbles  the  gaseous 
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phase  consists  almost  entirely  of  the  vapor  of  the  surrounding  fluid. 

Gas  bubbles  have  been  extensively  investigated,  particularly  in 
the  areas  of  surface  oscillations  and  mass  diffusion  effects.  However, 
for  the  field  of  cavitation  they  are  not  relevant. 

The  category  of  vapor  bubbles  is  subdivided  into  two  other  topics 
based  upon  the  extent  to  which  thermal  effects  alter  the  internal  pres¬ 
sure,  plt  and  thus,  bubble  behavior.  The  strong  pressure  dependences 
of  equilibrium  bubble  pressure  and  vapor  density  act  to  significantly 
reduce  p1  for  a  growing  bubble  as  evaporation  at  the  bubble  surface 
cools  the  interior.  Using  an  energy  balance,  and  assuming  that  heat  is 
supplied  by  a  liquid  layer  which  has  thickness  comparable  to  the  diffu- 
slon  length  (Dt)  ,  the  "thin  thermal  boundary  layer"  assumption,  this 
temperature  difference  can  be  estimated  as 

*P*(T)L 

AT  *  -  ~  -r -  (2) 

3(Dt)^  pc 

where  L  is  the  latent  heat  of  vaporisation,  p®(T)  is  the  equilibrium 
vapor  density  at  temperature  T,  D  is  the  thermal  dlffuslvity,  c  is  the 
heat  capacity  and  t  is  the  time  required  to  grow  to  radius  R.  For 
water  at  15°C,  with  R  ■  0.1  cm  and  t  ■  10"*  sec,  AT  ■  0.2#C,  while  at 
100°C,  AT  -  13°C.  For  the  former  case,  in  which  the  thermal  effect  is 
not  expected  to  be  Important,  the  proper  term  is  oaoitation  bubble.  In 
the  latter,  and  all  cases  where  thermal  effects  dominate  inertial  effects, 
the  result  is  'boiling  or  vapor  bubbles. 

The  cavitation  bubble  is  more  relevant  to  this  work  and  fortunately, 
is  the  simpler  case.  For  constant  p^»  and  neglecting  viscous  effects 
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where  py  -  1b  the  vapor  pressure.  pm  is  ambient  pressure  and  initial 
values  are  denoted  by  the  subscript  sero.  For  Rq  ■  0  and  neglecting  the 
surface  tension  term,  the  Rayleigh  result  for  time  to  complete  collapse 
Is 


Quantitative  experimental  confirmation  of  this  value  was  given  by 
Lauterborn  (1972).  He  was  able  to  generate  "empty"  cavities  by  focus¬ 
sing  a  pulse  from  a  Q- switched  ruby  laser  on  a  point  on  the  Interior  of 
a  liquid  mass.  Agreement  between  experiment  and  equation  (4)  was  ex¬ 
cellent,  despite  the  compressibility  of  the  real  liquid.  Note  that  this 
compressibility  Is  expected  to  become  significant  in  the  final  stages  of 
collapse,  when  the  Mach  number  of  the  bubble  becomes  large. 

The  large  magnitude  of  the  thermal  effect  on  boiling  bubbles  la 
Illustrated  in  Figure  2.  Here  the  theoretically  predicted  radius  vs. 
time  profiles  for  the  Plesset  analysis  is  shown  along  with  the  Rayleigh 
result  for  water  at  103*C.  Excellent  experimental  agreement  with  the 
Flasset  prediction  was  found  by  Dergarebedlan  (1953).  Ha  observed 


Figure  2. 

Thermal  Effects 
on  Spherical  Bubble 
Growth 

(Dergarabedi&n ,  1953 
Plesset  &  Zwick, 
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spontaneous  nucleation  in  superheated  water. 

Zwlck  and  Plesset  (1955)  have  also  analysed  the  collapse  of  a 
bubble.  The  distinction  between  cavitation  and  boiling  bubbles  is 
unimportant  here.  Despite  the  heating  which  occurs  with  vapor  conden¬ 
sation,  Figure  3  shows  little  deviation  from  the  Bayleigh  result. 

Flynn  (1975,  1976)  has  developed  a  more  complicated  model  includ¬ 
ing  compressibility  and  motions  within  the  cavity  as  well  as  heat  con¬ 
duction  and  viscosity.  His  results  cannot  be  summarized  easily,  but 
result  in  predictions  of  numerous  regimes  of  bubble  dynamics  which  have 
yet  to  be  experimentally  verified. 

The  results  most  relevant  to  the  present  study  suggest  that  the 
growth  and  collapse  of  a  spherical  cavity  in  water  are  usually  domi¬ 
nated  by  inertia  and/or  surface  tension,  not  by  viscosity.  The  small 
amount  of  viscoelasticity  expected  for  dilute  polymer  solutions  falls 
in  the  same  category  as  viscosity.  Both  viscosity  and  viscoelasticity 
berime  more  important  for  small  radial  oscillations  (Tanas swa  and  Tang, 
1970,  Zana  and  Leal,  1975,  Tang  and  Teh*  1983)*  but  are  still  minor. 

Experimental  testing  of  these  results  has  been  conducted  almost 
exclusively  on  aqueous  systems.  Spark  gap  and  laser  techniques  have 
been  used  to  concentrate  the  energy  necessary  for  cavitation  inception. 
In  general*  good  agreement  with  theory  has  been  found  (see  Figure  2). 


More  detail  can  be  found  in  the  section  entitled  "Experimental  Program". 

These  analyses  have  all  been  carried  out  under  the  assumption  that 
spherical  symmetry  of  the  bubble  is  maintained.  Plesset  (1956)  per¬ 
formed  a  linear  stability  analysis  on  the  nearly  spherical  interface 
between  two  immiscible,  incompressible,  inviscid  fluids.  Be  assumed  a 
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drop  shape  given  by 

'.urtac.  »•'>  '  *<'>  +  l  V*>V*>  (5> 

II 

where  is  a  spherical  harmonic  of  degree  n  (MacRobert >  1948; Hobson, 1955) , 

and  a  is  initially  saall  compared  to  the  equivalent  radius  R(t).  The  analy- 
n 

sis  proceeded  conventionally  in  order  to  determine  the  conditions  for  which 
aQ(t)  will  grow,  suggesting  the  shape  is  unstable,  and  those  for  which 
it  will  not,  implying  stability.  The  stability  condition  which  results 
can  be  expressed  in  terms  of  the  function  a  (t). defined  by 

U 

R  3/2 

*  (t)  -  (r^)  a  <t>  (6) 

n  K  n 

which  was  shown  to  be  governed  by 

#  R  “  G(t)a  -  0  (7 A) 

n  n 


where 


G(t) 


3  Ju2  R  r  3  n(n-l)pa-(n+l)(n+2)pi1 

2  V  +  R  L  2  +  np2  +  (n+i)pi  J 


(n-l)n(n+l)(n+2)o 
Inp  a  ♦  (n+2  )  p  i JR* 


(7B) 


and  o  is  the  surface  tension,  and  pi  and  pj  are  the  fluid  densities  of 
the  interior  and  exterior  fluids,  respectively.  The  details  of  the 
amplitude  profile  depend  on  the  particular  initiation,  but  a  general 
statement  can  be  made 


Si 
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G(t)  <  0  pronotes  stability 

G(t)  >  0  allows  instability 

this  means  that  surface  tension  always  has  a  stabilising  effect*  and 
from  Equation  (6).  increasing  R  or  bubble  growth  is  seen  to  have  a 
stabilizing  effect*  while  bubble  collapse  promotes  instability. 

Flesset  and  Mitchell  (1956)  then  performed  a  more  involved  analy¬ 
sis  for  a  vapor  cavity,  neglecting  the  density  of  the  vapor  and  the 
viscosity  of  the  vapor  and  the  liquid.  Their  results  showed  that  an 
expanding  vapor  cavity  is  stable*  i.e.  if  |an(0)|/RQ  «  1,  then 
|an(t)|/R(t)  «  1.  For  a  oollapeing  cavity*  distortion  amplitudes  re¬ 
mained  small  as  long  as  1.0  >  R/R  >  0.2*  but  as  R  0,  a  (t)  increased 

—  o  —  n 

in  magnitude  as  R  .  Thus,  the  spherical  shape  is  unstable  for  the 
later  stages  of  bubble  collapse.  These  trends  are  valid  even  when  a 
small  viscous  effect  is  included  in  the  treatment  (Prosperettl  & 
Seminars,  1978).  Viscosity  does  tend  to  damp  the  growth  in  amplitude 
of  the  higher  order  harmonics. 

To  determine  the  linear  stability  of  the  bubble  shape  only  the 
long  time  behavior  of  the  distortion  is  necessary.  The  bubbles  were 
assumed  to  exist  in  an  infinite  fluid*  quiescent  apart  from  the  effects 
of  the  dynamics  of  the  cavity;  this  is  a  spherically  symmetric  geometry 
giving  an  equilibrium  shape  of  similar  character.  Of  course,  this  is 
not  the  only  environment  in  which  a  bubble  may  arise.  There  may  be 
conditions  Imposed  on  the  fluid  away  from  the  bubble*  such  as  an 
elongatlonal  flow  or  a  solid  boundary*  or  the  bubble  may  be  initial¬ 
ised  as  nonspherlcal  and  a  detailed  description  of  its  evolution  de¬ 
sired.  Bubbles  have  provoked  some  study  in  this  regard,  but  more 
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often  liquid  drops  and  solid  bodies  have  been  examined,  giving  rise  to 
general  techniques  for  nonspherlcal  shapes. 

With  bubbles  and  any  other  initially  spherical  bodies  in  an  im¬ 
posed  flow,  the  primary  method  of  study  has  been  expansion  of  the 
velocity  profile  in  terms  of  spherical  harmonics  (Cox,  1969,  Shlma  & 
Nakajima,  1977).  Use  of  this  technique  is  usually  predicated  upon  an 
assumption  of  creeping  flow,  at  least  in  the  Immediate  vicinity  of  the 
bubble.  Under  this  restriction,  Happel  and  Brenner  (1965,  section  3.2) 
offer  a  practical  presentation  of  the  use  of  spherical  harmonics  in 
Lamb's  general  solution  (Lamb,  1945). 

The  presence  of  any  interface  in  the  vicinity  of  a  bubble  alters 
the  flow  induced  by  the  cavity  dynamics  from  a  spherically  symmetric 
velocity  field.  For  a  solid  wall  the  collapse  process  results  in  a  jet 
impinging  on  the  wall  which  may  be  one  of  the  major  causes  of  the  de¬ 
structive  action  of  cavitation.  Experimentally  the  jet  has  been  ob¬ 
served  in  laser-induced  bubbles  (Lauterborn  &  Bolle,  1975)  in  qualita¬ 
tive  agreement  with  numerical  simulations  (Plesset  &  Chapman,  1971). 
(See  Figure  4.)  This  type  of  behavior  can  be  predicted  using  either 
integral  methods,  employing  the  Bernoulli  equation  while  neglecting 
viscous  stresses  (Voinov  &  Voinov,  1975,  1976),  or  marker-and-cell 
simulations  (Mitchell  &  Hammitt,  1973),  which  can  include  viscosity. 

The  viscous  effects  for  water  are  negligible,  dominated  by  surface 
tension  and  Inertia.  (See  Table  2  for  representative  values.)  The 
most  Important  parameter  appears  to  be  the  distance  of  the  bubble  from 


the  wall 
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Order 


TABLE  2 

of  Magnitude  of  Terms  In  Radius  Equation 
Inertial,  Surface  Tension  and  Viscous 

R  3  x  10-1  cm 

p  1  gm/cm3 

V  10“2  gm /  (cm-sec) 

o  'v*  60  gm/sec2 

At  *  10“ 3  sec 

«=>  R  ~  R/At  ^  3  *  102  cm/sec 

Inertial:  p  R2  10s 

2o 

Surface  Tension:  ^  400 

Viscous:  y  —  %  10 
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When  the  Interface  is  not  a  solid  vail,  hut  is  defornable  like  a 
free  surface,  an  "entirely  new  form"  (Gibson  &  Blake,  1981)  of  collapse 
can  be  observed  experimentally.  Depending  on  the  physical  properties 
of  the  boundary  —  'inertia*  and  'stiffness*  —  the  bubble  nay  nigrate 
to  or  from  the  interface  as  the  bubble  collapses  or  pulsates.  A  jet 
forms  and  moves  in  the  same  direction  as  the  overall  siigratlon. 

An  initial  nonsphericity  is  important  even  in  the  absence  of  an 
imposed  flow  or  solid  boundary.  The  requisite  mathematics  has  been  put 
forth  by  Hsleh  (1965,  1974),  who  includes  equations  for  heat  and  mass 
transfer  effects,  compressibility  and  also  variational  methods.  Sim¬ 
plified  numerical  results  by  other  authors  (Chapman  &  Plesset,  1971) 
show  that  an  Initially  prolate  spheroid  (Figure  5A)  may  form  -wo  jets 
upon  collapse,  while  an  oblate  spheroid  (Figure  5B)  leads  to  a  dumbbell 
form.  The  most  relevant  feature  of  these  nonlinear  results  is  that  the 
linear  theory  of  Plesset  and  Mitchell  (1956)  agrees  well  until  the  final 
stages  of  collapse. 

Very  little  has  been  offered  here  so  far  about  non-Newtonian  or 
viscoelastic  effects,  either  generally  or  in  relation  to  bubbles.  In 
the  general  view  it  is  best  to  refer  to  the  numerous  treatises  on  the 
subject  (Bird,  Armstrong  &  Hassager,  1977,  Schowalter,  1978).  The  most 
important  general  consideration  is  that  fluids  which  exhibit  viscoelastic 
effects  may  behave  in  a  manner  which  runs  counter  to  an  observer's 
"Newtonian"  intuition  (Schowalter,  1978),  e.g.  the  Welssenberg  rod¬ 
climbing  effect.  Furthermore,  the  history  or  pre-treatment  of  such  a 
material  can  be  vitally  Important  to  its  performance.  In  modelling 
there  ere  innumerable  complications  which  almost  invariably  force  a 


compromise  between  rigorously  formulated  constitutive  equations  of  fluid 
behavior  and  simpler  forms  which  result  in  tractable  mathematics.  The 
applicable  models  and  quantitative  results  are  also  highly  dependent 
on  the  particular  flow  present*  although  qualitative  trends  are  usually 
similar.  These  difficulties  must  be  considered  in  the  quasl-llnear 
model  employed  in  the  present  work;  careful  attention  is  given  to  pos¬ 
sible  nonlinear  effects.  (See  Appendix  B) 

Dilute  solution  non-Newtonian  effects  on  bubbles  brings  attention 
back  to  cavitation  inhibition.  Spherical  bubbles  show  little  effect 
(Ting,  1975,  1977,  Ting  &  Ellis,  1974),  experimentally  or  theoretically. 
Experiments  involving  nonspherlcal  bubbles,  by  Chahine  (1981)  and  others 
(Chahlne  &  Fruman,  1977,  Gibson  &  Blake,  1981)  with  collapsing  bubbles 
near  interfaces,  show  a  "significant"  delay  in  the  creation  of  the 
micro jet  when  250  ppm  of  polyox  are  added  to  the  solution  and  the  inter¬ 
face  is  solid.  Observations  near  a  free  surface  also  show  an  effect 
from  dilute  additives.  However,  these  results  are  not  definitive  and 
the  relationship  between  them  and  cavitation  effects  is  still  a  matter 
of  speculation.  Observations  by  different  workers  may  even  appear  con¬ 
tradictory.  Chahlne  suggests  that  the  level  of  nonsphericity  in  bubbles 
is  reduced  upon  addition  of  polymer  solute,  while  Oba  (1978)  reports  the 
opposite  effect  in  his  flow  experiments. 


III.  THE  EQUATIONS  OF  MOTION 


1.  Spherical  Dynamics 

For  the  cavitation  model  developed  here  the  nonspherlcal  shape 
will  be  coupled  with  the  overall  growth  or  collapse  of  the  bubble. 
This  volume  change  is  treated  through  an  analysis  of  the  equivalent 
spherical  bubble.  Continuity  requires  that  the  velocity  field  in  an 
incompressible  fluid  in  this  symmetrical  situation  have  the  form 


^spherical  (r,t^r 


Rjt), 

r2 


(1) 


in  spherical  coordinates  with  origin  at  the  bubble  center.  R(t)  and 
R(t)  are  the  Instantaneous  radius  and  its  time  derivative,  respectively. 
The  fluid  is  assumed  to  be  quiescent  at  infinity  so  the  angular  velocity 
components  are  Identically  zero. 

The  r-component  equation  of  motion  integrated  from  the  surface 
r  ■  R  to  infinity  then  gives 


p[RR  +  |  R2]  -  P£-Pa  + 


(2) 


where  p,  the  liquid  density,  t,  the  extra  stress  tensor,  P£,  the  liquid 
pressure  at  r  •  R,  end  Pg,  the  ambient  pressure, must  now  be  specified. 
By  neglecting  the  viscosity  and  density  of  the  internal  ges  the  bubble 
is  assumed  to  have  a  uniform  internal  pressure.  This  should  be  accu¬ 
rate  as  long  as  the  ratio  betwesn  internal  end  liquid  viscosities  is 
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small  and  the  Internal  Mach  number 


^internal  '  »/(3  *  10*  «  1 


(3) 


since  the  speed  of  sound  in  air  is  approximately  equal  to  this  denomi¬ 
nator.  Surface  tension  can  be  Included  so  that  Equation  (2)  can  be 
expressed  as  a  normal  force  balance  across  the  surface 


p[RR  +  |  R2]  -  Pg(t>  +  Py(t)  -  Pa(t) 


(4) 


where  P  is  the  Internal  pressure  of  non-condensable  gas.  Py  is  the 
© 

vapor  pressure  of  the  liquid  at  the  surface  temperature  and  o  is  sur¬ 
face  tension.  This  form  of  the  equation  is  sufficiently  general  to  be 
capable  of  including: 

A)  The  nature  of  the  non-condensable  gas 

e.g.  Pg(t)  -  Pg(R)  -  Rg  (-|)3t  (5A) 

for  a  polytropic  gas 

B)  Thermal  effects 


Pv(t)  -  Py(T),  o(t)  -  o(T) 
and  v(t)  -  v(T) 


(5B) 


C)  Different  constitutive  stress  expressions 
e.g.  rrr  •  f(Ix.  12.  I»)  for 
strain  invariants  I 


(5C) 


The  simplest,  useful,  non-taro,  stress  expression  is  that  for  a  purely 


Newtonian  fluid  with  constant  viscosity  p.  For  this  velocity  field 
the  radial  deformation  rate  Is 


rr 


(r,t)  - 


«) 


so  that 

R 

Non-Newtonian  expressions  may  be  complicated  by  higher  time  deriv¬ 
atives  or  integration  over  past  times  combined  with  material  objectivity 
constraints.  When  the  time  Integration  is  performed  in  a  Lagrangian 
frame,  following  each  fluid  element,  objectivity  is  satisfied  for  this 
simple,  symmetrical  flow  situation.  The  Lagrangian  coordinates  can 
easily  be  related  to  the  Eulerian,  laboratory-fixed,  coordinates  by 

(r')s  -  rs  +  R3(t’)  -  R*(t)  (8) 

where  r'  denotes  the  position,  at  the  past  time  t*»  of  a  particle  which, 
at  present  time  t,  resides  at  radius  r.  Time  integration  must  be  per¬ 
formed  at  constant  r*. 

A  mathematically  simple  model  which  includes  stress  accumulation 
with  fading  memory  was  employed  by  Fogler  and  Goddard  (1970,  1971),  who 
specified  a  relaxation  modulus  (memory  function)  N(t)  such  that 

I 
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Trr(t)  "  {  Mt-Oe^t^dt* 

— CO 

N(t)  -  p  act)  +  Ge  exp  (-t/X) 


(9) 


when  6 (t)  is  the  Dirac  delta  function.  The  elastic  contribution  is 
characterised  by  two  parameters:  X,  a  relaxation  time*  and  Gq,  an 
elastic  modulus.  In  the  limit  of  xero  elasticity  a  Newtonian  fluid 
with  viscosity  y  is  described.  In  the  Lagranglan  frame*  the  stress 
term  becomes 


GO 


R 


t 

-  -4  j  N(t-t’) 


R(t,)R2(t,)i^(|^p-)dt, 
R*(t')  -  R* (t) 


(10) 


A  similar  result  obtains  for  an  Oldroyd  three-constant  model. 
Ting  (1977)  shows  that  the  equivalent  elastic  modulus  is  expected  to 
be  directly  proportional  to  polymer  concentration.  To  generate  a  rad¬ 
ius  profile  initial  conditions  need  to  be  specified.  Ting  chose  to 
Impose  equilibrium  before  time  t  •  0,  then  postulated  a  step  change  in 

ambient  pressure  of  magnitude  P*  at  time  t  •  0  to  initiate  notion.  By 

* 

his  convention,  when  P  <0  the  bubble  grows.  The  numerical  solution 
predicts  changes  in  R(t)  values  of  less  than  2Z  between  pure  water  end 
a  slightly  viscoelastic  medium. 

Fogler  and  Goddard  present  lsrge  elastic  effects*  i.s.  changes 
in  the  R(t)  profiles,  but  for  psrsmeter  values  which  minimise  surface 


„rt»> 


tension  and  correspond  to  Gq  such  beyond  the  snail  values  expected  by 
Ting.  Since  the  present  work  was  activated  by  dilute  solution  phe¬ 
nomena,  values  closer  to  Ting's  were  investigated.  This  neant  little 
alteration  of  the  steady  shear  viscosity  was  allowed.  One  wanner  by 
which  to  isipose  this  requlreaent  was  to  restrict  the  relative  values 

of  G  and  X  by 
o 


where 


GoX  -  v 


v  *  0  <10”s  Pa-s) 


which  gives  an  elastic  contribution  to  the  steady  shear  stress  equal  to 
the  viscous  component.  In  calculations  the  initial  physical  parameter 
values  were  chosen  to  correspond  to  superheated  water  at  103*C,  since 
experimental  data  on  spherical  bubble  dynamics  were  available 
(Dergarabedlan,  1953). 


TABLE  3 

Parameter  Values  for  Initial  Calculations 


Pv  -  845  n  Hg 
Pa  -  760  an  Hg 


o  ■  58.9  dyne  cm 
p  -  0.961  gm/cm* 


-10*  dyne/cm2  <  P*  <  10s  dyne/cm2 


10”*  sec  <  X  <  1.0  sec 
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For  all  P*  and  X  values  the  radius  profiles  changed  by  about  IX  or  less 

between  a  purely  viscous  and  viscoelastic  fluid  subject  to  condition  (11). 

* 

As  the  magnitude  of  P  increased,  inertia  became  more  important  since 
velocities  and  accelerations  became  larger.  For  more  detailed  results, 
see  Section  VI. 


2.  Viscous  Nonsphericities 

The  relaxation  of  the  assumption  of  spherical  symnetry  greatly 
complicates  the  analysis,  even  for  the  purely  viscous  case,  since  0  and 
<J>  components  and  dependencies  now  exist.  A  general  analysis  might  begin 
with  the  full  equations  of  motion  in  spherical  coordinates.  However, 
the  Navler-Stokes  equations  include  nonlinearities  from  the  convected 
inertial  terms.  These  introduce  computational  difficulties  which 
are  compounded  for  viscoelastic  fluids  by  nonlinear  co-deformatlonal  or 
co-rotational  constitutive  models. 

Prosperettl  (1977 ,  1980)  neglected  the  non-linearities  in  the 
purely  viscous  case  by  imposing  the  condition  that  all  departures  from 
spherical  symmetry  be  small.  The  shape  of  the  bubble  surface  J( r,0  4>,t) 
can  then  be  expanded  in  spherical  harmonics. 


I  (r,6,*,t): 

r  -  R(r)  -el  a  (t)V®  (0,*)  -  0 

it  u 


(12) 


for  some  small  parameter  e  which  must  be  determined  to  allow  the  desired 

mathematical,  and  physical,  simplifications.  Here,  a  (t)  is  the 

n 


magnitude  of  the  spherical  harmonic  ^(6.*).  The  spherical  harmonice 
determine  a  particular  angular  dependence  and  are  defined  as  the  aolu 


tions  to 


[n(n+l)  +  |q  sin®  g0  +  iiiFe  “  0 


or,  equivalently 


(e »«p)  -  ei|m'*  p (cose) 


(13) 


|m  f  v 

for  the  Legendre  polynomials  Pfl  (x) . 

Of  course,  the  full  velocity  field  V(r,6.*,t)  must  satisfy: 


H  .  |I  +  (V  .  V)V  -  ^  div.  0 
Dt  3t  ~  P 


(14) 


where  the  stress  tensor  a  for  a  Newtonian  fluid  in  Cartesian  coordinates 
is 

/su  av  \ 

(15) 


m  9VA 
’ij-  '  *  eU  *  V\»T  + TtJ 


The  surface  .hap.  [.  ebove,  is  con.lst.ut  «itb  .  perturbation  expansion 
of  V  and  the  pressure  field  P  as 


y(r,e,*.t)  -  ▼  +  t  v  +  t  y  +  0(c2) 


(16) 


P(r,0,4,t)  -  P0  4  t  pp  +  e  Pv  +  0(e2) 


(17) 


when  v  and  p  obtain  for  the  unperturbed  spherically  sy»etric  case 


previously  presented 


PQ  -  Pa(t)  +  p[(R2K  +  2RR2)/r  -  |  (18) 

The  first  order  asymmetries  can  be  divided  into  two  parts;  the 
first  arises  in  an  inviscid  potential  analysis.  Plesset  (1954)  shows 

v  «  V  6  (19) 

-P  P 

o+2  • 

S  -  -  sr  Tfi Il(t)  +  2*«> 

and 

Pp|  (R(t)a(t)  +  3R(t)a(t)  +  2R(t)a(t)]T^(e,*)  (20) 

r-R 

The  second  type  of  0(e)  term  is  necessitated  by  the  introduction 
of  viscous  dissipation  and/or  elastic  energy  storage.  Since  potential 
functions  such  as  the  inviscid  velocity  profile  have  sero  curl,  this 
second  velocity  can  be  treated  through  the  vorticity  w 

c  w  Jxy  -  c  Jxyv  (21) 

which  is  governed  by  the  usual  vorticity  equation  for  the  purely  viscous 
case.  At  first  ordsr  in  c 

♦ 

du 

j~  +  Jx(«  *  yD)  •  -v  V  *(v  x  w)  (22) 


for  kineaatlc  viscosity  v 
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The  usefulness  of  this  analysis  in  terms  of  vortlcity  is  not  only 
the  elimination  of  the  need  to  explicitly  determine  the  pressure  field, 
but  also  a  decomposition  of  $  is  possible  since  it  is  solenoidal 
(V  *  u  -  0).  The  vortlcity  can  be  expressed  as  the  sum  of  two  other 
vector  fields,  one  poloidal  and  the  other  toroidal,  which  can,  by  defi¬ 
nition,  be  expressed  through  two  series  of  eoalar  functions. 

-  §  +  T  (23) 

S(r,e,*,t)  -  l  l  2xVx[Sn(r,t)YjJ(e,*)§rl  (24) 

n  m 

T(r,e,*,t)  mH  ?x[Tn(r,t)Y^(e,*)er]  (25) 

n  m 

Since  these  two  vector  fields  are  orthogonal,  the  linearized 

equation  can  now  be  written  as  two  series  of  independent  scalar  equations. 

For  each  index  n: 

f +R(f)|f-  vt0-n(n+1>f2]  (26) 

f  +  h  -  »«"*>  7*1  <27> 

The  properties  of  poloidal  and  toroidal  fields  (Chandrasekhar, 


App.  Ill)  can  next  be  employed  to  show  that: 


(28) 


*,  ■  I  +  -  V*,  \ 

n,m  * 

Continuity  considerations  Tequire  the  presence  of  the  final  tern  which 
mi«t  satisfy 

V2*  -  7  •  T  Ve  (29 A) 

n  n  tl~x 

Integrating  this  equation 

♦n  '  ^  KCt)  +  551  f 

R 

+  I^r  *2n+1«»(t>  +  ssr  f  s"tlTn<s.t)ds]r-(,,+1)}  <29B) 

R 

where  the  "constant"  of  integration  a  (t)  is  determined  by  the  ambient 

n 

flow  conditions.  When  the  flow  decays  at  infinity, 

OB 

«„<»  •  -  iSr  I  s'“  Vs-t)ds  •  <*» 

R 

These  defining  equations  for  S  and  T  and  their  relationship  to 
the  perturbation  velocity  y^  allows  an  a  posteriori  evaluation  of  e. 
The  simplification  of  these  equations  is  only  valid  when  the  nonlinear 
terms,  corresponding  to  convective  terms  in  the  substantial  time 
derivative,  which  are  6(e2),  are  much  smaller  than  the  6(c)  terms. 
These  restrictions  are  considered  in  Appendix  B. 
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Boundary  conditions  on  S  and  £  are  now  required  to  specify  the 
solution  (Cox,  1969,  Scriven,  1960).  For  an  invlscld  and  massless 
internal  gas  the  velocity  profile  oust  satisfy  the  zero  tangential 
stress  condition  and  the  balance  of  normal  stresses  at  the  interface. 
The  former  condition  becomes: 


|^-  {r“2[yS(r,t)]}r-R  -  0  (31) 

2  (2n+l)Rn~2  [-pa(t)  ]  +  -^~^[pT(R,t)] 

•  e 

-  2(n+2)w  ~  -  (n-l)pa  (32) 


where  the  amplitude  a(t)  and  its  first  derivative  enter  the  equation 
through  considerations  of  surface  curvature.  The  viscous  contribution 
to  pressure  in  the  external  fluid  at  r  -  R  is 


-  npY®  { 
r-R  *■ 


v  + 
R 


(*) 

V 


Rl(f)9-l)(f)nT(S,t)ds|  (33) 


and  it  can  be  used  to  express  the  normal  stress  boundary  condition  at 
r  -  R. 


+  l^+r  * '  I  <*«)CB-i>v]i(t) 


+  l-  R  +  2(n-l)(n+2)v  |2  +  (n-l)(n+2)  ^]a(t) 
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■i 

% 

n 

\ 

i 

'1 

( 


+  n(n+2)v  -  n(|)  J  [1  -  <|)*]  (|)nT(S,t)dS  (34) 

R 

Since  the  bubble  shape  determined  by  the  a  (t)  is  of  primary  interest 

U 

this  last  equation  may  be  viewed  as  the  principal  expression  of  interest. 
It  is  an  equation  of  motion  for  bubble  shape,  while  the  other  two  are 
boundary  and  consistency  conditions. 

Prosperetti  solved  only  for  the  case  of  constant  bubble  radius. 
Since  the  coefficients  of  a(t)  and  its  derivatives  are  then  constants 
and  only  the  inhomogeneous,  forcing  function  is  time-dependent,  solution 
of  the  system  is  considerably  simplified.  Laplace  transform  techniques 
can  be  applied  to  determine  the  asymptotic  behavior  of  a^Ct)  analyti¬ 
cally.  However,  to  generate  a  complete  an(t)  profile  numerical  tech¬ 
niques  are  necessary  to  invert  the  solution  in  the  Laplace  domain.  Of 

course,  initial  conditions  such  as  the  initial  amplitude  a  (0)  and  non- 

n 

spherical  velocity  an(0)  are  also  necessary. 

Results  for  this  particular  case,  where  there  is  no  motion  pre¬ 
vious  to  time  t  ■  0,  show  nonspherlcal  amplitudes  undergoing  damped 
oscillations.  Initially,  the  frequency  is  determined  by  surface  ten¬ 
sion  and  damping  factor  by  viscosity,  but  no  such  simple  statements 
can  be  made  about  the  behavior  at  long  times. 

These  Laplace  techniques  are  not  applicable  when  the  equivalent 
radius  R  is  not  constant,  but  the  results  for  this  case  will  probably 
be  more  interesting  in  view  of  the  stability  analyses  summarised  in 


i 


Section  II  and 


"Indeed,  it  ie  expected  that  in  a  flow  with  ao 
strongly  converging  (or  diverging)  streamlines 
as  the  spherical  one  the  accompanying  concentra¬ 
tion  (or  dilatation)  of  vorticity  would  play  an 
important  role  in  the  process  of  energy  dissipa¬ 
tion  and  hence  in  the  equation  of  motion  for  a(t)w. 
(Prosperetti ,  1977,  pp.  345,346) 

Prosperetti  (1980)  has  also  approximated  the  solution  of  some 
cases  of  variable  R  by  making  a  series  of  approximations.  He  concludes 
that  small  viscosity  will  not  substantially  inhibit  the  instability  of 
the  spherical  shape  for  moderate  n.  However,  it  may  change  the  mode 
of  breakup  if  the  bubble  collapses  into  many  fragments. 

Before  proceeding,  the  neglect  of  the  pololdal  field  $  should  be 
justified.  All  equations  governing  this  function  (Equations  26  &  31) 
are  homogeneous,  so  it  is  never  present  unless  initiation  requires  its 
existence.  Physically,  a  zero  pololdal  field  corresponds  to  an  absence 
of  circulation  or  bulk  rotation  (Batchelor,  1970  p.  93)  in  the  fluid. 
Mathematically,  pololdal  fields  are  necessary  when  the  velocity  con¬ 
tains  components  with  angular  dependences  given  by 

V  x  [iJJ  (6,4)  gr]. 


3.  Viscoelastic  Nonsphericities 

Extension  of  this  perturbation  technique  to  include  viscoelas¬ 
ticity  is  straightforward  in  principle,  however  the  level  of  complexity 
of  the  resulting  equations  is  much  higher  and  there  are  additional 
ramifications  of  the  linearisation.  An  integral  model  (Equations  9  6  11) 
of  fluid  behavior  has  been  chosen  in  part  because  explicit  expressions 
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for  c't.'c  stresses  are  Chen  available.  This  introduces  an  integro- 
dlfferential  form  to  the  system.  Material  objectivity  requirements , 
e.g.  co-deformatlonal  or  co-rotational  integrations,  will,  in  general, 
introduce  nonlinearities  Which  must  be  evaluated  (See  Appendix  B.2). 

The  extra  stress  expressions  in  the  purely  viscous,  Newtonian, 
case  have  the  relatively  simple  form  of  the  terms  in  equations  (27), 
(32)  and  (34)  which  Include  a  viscosity  constant  y  or  v.  Any  further 
complications  are  removed  by  the  instantaneous  nature  of  the  process 
of  viscous  dissipation  and  its  characterisation  here  by  a  single,  time- 
and  position-independent  physical  constant.  Any  spatial  variation  in 
the  viscosity  would  disallow  the  amplifications  which  result  from  the 
manipulation 


V  •  y[(V?)  +  (VV)T]  -  2yV2V  (35) 

and  others  like  it.  To  include  viscoelasticity,  all  these  extra  stress 
terms  must  be  altered  with  careful  consideration  of  the  order  of  opera¬ 
tion  and  reference  frame  or  coordinate  system  in  which  they  are  performed. 

The  fundamental  difference  when  viscoelastic  stresses  are  calcu¬ 
lated  is  the  need  to  operate  not  merely  on  the  instantaneous  strain 
rate,  but  also  on  its  appropriate  time  derivative  in  a  differential 
model  or  on  the  pertinent  past  values  for  an  integral  equation.  Follow¬ 
ing  Fogler  and  Goddard, a  Maxwell-type  Integral  model  is  employed  here 
since  it  is  tenable  with  superposition  of  simple  viscous  and  fading 
elastic  effects  and  la  the  small-deformation  limit  of  many  other,  more 
general  expressions.  Beginning  with  velocity  field  (28)  and  boundary 
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conditions  (32)  and  (34)*  terns  of  the  form 


v  j  (R;r,t) 


(36) 


and  replaced  in  the  linear  theory  by 


J  N(t-t')  y  (R* ;r' ,t')dt' 


(37) 


This  is  not  a  simple  Bvlerian  Integral  over  past  times.  Rigorous  mater¬ 
ial  objectivity  constraints  would  require  a  Lagranglan  frame  with  co- 
rotatlonal  or  co-deformatlonal  contributions  to  take  convective  and 
orientational  changes  in  inter-fluid  element  spatial  configurations 
into  account. 

Radial  convection  is  incorporated  to  0(e°)  just  as  it  was  in  the 
spherical  case.  A  new  radial  position  coordinate  H  is  defined  here 
which  is  equivalent  to  the  Langranglan  coordinate  (8)  for  symmetric 
dynamics.  The  variable  transformations  are  defined  by 


H 


.*  _  »3 


(38) 


T  -  t 


so  that 


T(H,t)  -  T(r.t) 


(39) 


Since  the  T-equatlon  results  directly  from  the  moment  us  or  vor- 
ticity  equations  the  right  hand  side  (IBS)  of  this  equation  must  be 
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derived  from 


(?*)? 


T 

•  |  K(t-T')Y(H, 


T’)dT' 


for  the  strain  rate  tensor  v.  The  new  T  equation  becomes: 


(40) 


|I  .  2K*i  5,  -  vL{2  ||  ♦  L‘  0  -n (n+1 )  |r> 

+  Go{i  f  «p(Ir1)1l(H,T,)<iT, 

—  00 

T 

+  2L2  J  exp(^L)I2(H,T,)dT’ 

»oo 

Ls  |  exp(^^,)Ij(H»T' )d't'  } 


(41) 


where 


It  -  2n(n-l)(n-2)Xn(H,T)L 


(n-2) 


-  2(n+l)(n+2)(xrt-3)Zn(H,t)L 

T  (H,t) 

+  n(n+l)  ~  ; - 


~(n+3) 


T  (H.O 

I2  -  -2(nz+n+4)  — — 


-  (2n-5)(n-l)(n-2)Xn(H,T)L 


(n-5) 


+  (n+2)(n+3)(2n+7)Zn(HtT)L 


-(.n+6) 
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T  (H*t) 

I3  -  2(n24tt+8)  — £7 - 


+  2(n-l)(n-2)(n-5)Xn(H,T)L 


(n-8) 


2 (n+2) (n+3) (n+8) Zn (H. t)L 


l!%  +  4  5i 

L  3  2  L*  3H 


-(n+9) 


L(H)  -  C3H  +  R*)1/3 


VH*T)  •  °n(T)  ■  An(T)  +  2»fr 


1  s'\ 


(H,x)dS 


z„<H>t)  ■  SCT  R2n+llan(T)  -  An<T)I 


L, 

♦sSrj  ^  In(H,T)dS 

R 

and  the  function  A(t)  ia  zero  for  an  externally  quieacent  fluid  and 
is  introduced  in  the  far-field  conditions  (Section  IV).  See  Appendix 
A  for  a  derivation  of  this  expression. 

There  ere  elso  nee  complications  in  the  formulation  of  the  boundary 
conditions  et  the  bubble  surf see.  A  distinction  emerges  between  terme 
which  erlee  from  the  stress  expression  end  those  in  the  unit  normal 


™hb» 
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vector.  In  the  viscoelastic  case,  terns  obtain  such  as 


T 


and 


R(t) 


T 


N(t-t') 


R(t') 

RCP) 


dr’ 


(43) 


which  are  equal  in  the  viscous  Unit,  but  must  be  distinguished  in  this 
application.  The  normal  stress  boundary  condition  is  altered  both  in 
this  regard  and  by  changes  in  the  T-equation.  The  0(e)  v-pressure  is 
calculated  by  the  spatial  integration  of  the  equation  of  motion  (see 
Appendix  A) 


Pv(H«0,t) 


1 

j 


N(t-t') 


(44) 


{‘  iff  +  3  I*  -Ztte-iX**-2)  +  (n+2)(n+3)(^r)] 


R(n"5)(o-A))dT’ 


+ 


|  N(t-t'){-  3(|) 


+  2[(n-l)(n-2)  +  (n+2) (n+3) ]R(°~2 ) («-A) )  dr’ 
+  2^«2  Rn(R)(o.A) 


+  R  T(H-0,x) 


-  R 


f  •- 


_n  f2n+U  3A 
"  R  *  n+l;  3t 


The  new  normal  stress  boundary  condition  becomes  (see  Appendix  A) 


Ra 

n+1 


+ 


3 

n+1 


R k  -  2(n-l) (n+2) 


|  N'fdx- 


+  (-  gl)Ra  -  4(n-l)  (n+2) 


K* 


aR 

W 


dr' 


+  6n(n+l)  ^ 


f  N’(|)dr’ 


-  12  aR2  |  N'^dx'  +  (n-l)(n+2)  ^  a 
— «0 

T 

+  n(n+l)  J  H’(f>dT'  +  O  -  0 


(45) 


o  !■  the  surface  tension  constant  and  S'  ■  H(t“T*)*  mhich  is 
▼lsvsd  ss  an  equation  of  notion  for  a (t ) •  This  suit  be  solved 


I 


consistent  with  the  T-squstion  snd  the  ten gent iel  stress  boundary  con¬ 
dition.  (See  Appendix  A.) 
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|  H*  {2(0+2)  {J  + 


-  (n+1)  |+  2(2n+l)R(n"2)(o-A)J  dr' 

T 

-  6 (n+1)  |  |  H*  |  dr'  -  0 


(46) 


It  is  useful  to  note  that  these  equations  do  reduce  to  their  viscous 
counterparts  when  the  elastic  contribution  is  sero.  In  particular,  the 
spatial  Integrals  of  the  T  field  are  seen  to  cancel  identically  in  the 
Instantaneous  dissipation  linit.  (See  Appendix  A. 3.) 

Viscoelastic  calculations  for  the  considerably  simpler  special 
case  when  R  -  Rq,  have  been  performed  by  Inge  and  Bark  (1981)  using 
Laplace  transform  techniques  while  including  a  linearized  3-constant 
Oldroyd  model.  Their  results  supply  s  useful  check  on  techniques  which 
allow  time-dependent  R(t)  values. 

The  system  of  equstions  governing  the  asysnetric  dynamics  de¬ 
veloped  hers  is  linear  in  all  of  the  0(c)  quantltlea,  l.a.  an(t), 

T  (H,t),  A  (t)  and  their  integrals  and  derivatives.  This  linearity 

u  u 

emerges  from  the  assumption  that  terms  which  are  aecond  order  in  c» 
representing  conveeted  moment**  fluxes  and  also  convected  and  co- 
rotational  and/or  co-deformetional  eontriubtlons  to  the  non-Hewtonian 


i 
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stress,  are  small  enough  to  be  neglected.  However,  the  constitutive 
model  employed  does  not  neglect  all  nonlinear  contributions  to  the 
stress,  the  Lagrangian  coordinate  tranaformation  incorporates  convec¬ 
tion  by  terms  of  zero  order  in  c.  Thus,  this  is  not  merely  a  linear 
viscoelastic  model.  Furthermore,  at  this  first  order  in  e,  the 
integral  model  of  Equation  (III. 9)  is  equivalent  to  many  others  in 
integral  form  including  the  linear  Jeffreys  model  (Bird,  Armstrong  & 
Hassager,  p.  279),  the  co-rotatlonal  Jeffreys  model  (p.  328)  and  the 
Oldroyd  8-constant  model  (p.  371). 

The  neglect  of  terms  of  second  order  in  c,  restricts  the  condi¬ 
tions  for  which  the  equations  are  valid  to  motions  exhibiting  small 
deformations  from  the  purely  radial,  0(e°)  motion.  The  prescription 
for  "small"  values,  which  determines  c,  can  also  be  viewed  in  terms  of 
the  harmonic  expansion  which  is  the  basis  for  the  separation  of  angular 
and  radius-temporal  dependences  in  equation  (25).  By  treating  each  set 
of  equations  (41),  (44),  (45)  and  (46)  for  each  value  of  the  index  n. 
Independent  of  other  n  values,  the  linear  Independence  of  the  modes 
(spherical  harmonic  angular  dependences)  is  Invoked.  This  implicitly 
means  than  any  asyunetric  field,  nonsphericity,  velocity,  pressure  or 
stress  with  specific  angular  dependence  T^(0,$)  will  "excite"  and  be 
"excited"  by  only  quantities  with  that  same  angular  dependence.  In 
reality  this  can  be  violated  by  Inertial  terms  or  by  convected  elastic 
contributions  from  angular  or  radial  strains  and  the  resulting  stresses. 
In  the  viscoelastic  stress  term  the  possible  contributions  can  be 
evaluated  mathematically  and  this  is  done  in  Appendix  B.2. 
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IV.  INITIATION  AND  FAR-FIELD  CONDITIONS 


In  order  to  solve  for  aQ  and  thus,  the  shape  of  the  cavity,  the 
viscous  equations  (111. 27),  (111.32)  and  (111.34)  or  the  viscoelastic 
system  (III. 41),  (III. 44),  (111.45)  and  HI. 46)  must  be  solved  simul¬ 
taneously  subject  to  an  6(e°)  R(t)  profile  such  as  one  generated 
through  equation  (III. 4).  Specification  can  also  be  made  of  ambient 
flow  conditions  In  conjunction  with  (III. 29)  and  (111.30).  The  New¬ 
tonian  system  is  mathematically  simpler  than  the  equations  containing 
the  Maxwell-type  model  and  complete  problem  specification  can  be 
achieved  by  two  physical  parameters  y/p  and  o/p  and  two  shape  values 

per  mode,  a  (0)  and  a  (0)  at  6(c).  For  a  bubble  in  a  fluid  with  no 
n  n 

externally  Imposed  flow  the  system  is  well-posed. 

Viscoelastic  fluids  add  a  non-instantaneous  component,  caused  by 
fluid  memory  and  represented  here  by  time  integrals,  which  give  elastic 
stress  accumulation.  Behavior  prior  to  some  initial  time,  e.g.  t  -  t 
■  0,  is  relevant  to  the  determination  of  subsequent  dynamics.  It  is 
no  longer  sufficient  to  specify  ®n(0)  and  8^(0)  only;  the  memory  inte¬ 
grals  must  also  be  initialised. 

Time  zero  is  chosen  to  be  that  instant  at  which  6(e°)  spherical 
dynamics  are  initiated  i.e. 


R(t)  -  R<t) 
R(t)  -  Rq 


(1) 


so  from  equation  (1X1.3) 


R 


o 


P°  -  P°  -  P° 
g  v  a 

2o 


(2) 


This  also  initiates  the  value  of  the  elastic  contribution  up  to  time 

aero  as  aero  at  0(c°).  For  this  vapor  cavity  (P°  *  P  (t)  -  0),  notion 

8  8 

is  begun  by  a  step  change  in  ambient  pressure 


P  (t)  «  P°  -  P*  H(t)  (3) 

&  S 

where  H(t)  is  the  Heaviside  step  function.  The  sign  of  P*  determines 
the  type  of  spherical  dynamics 


p* 

<  0 

R  <  0 

Collapse 

* 

• 

p 

-  0 

R  -  0 

Static 

p* 

>  0 

R  >  0 

Growth 

For  a  particular  value  of  P* ,  and  the  physical  parameters  of  Table  3 
(in  Section  II),  R(t)  vs.  t  profiles  can  be  generated  from  equations 
(III. 3)  and  (III. 9)  with  (III. 11)  and  one  final  condition 


R(0+)  -  0 


(5) 


The  initial  conditions  on  the  nonspherlcal  terms  which  are  analo¬ 


gous  to  those  on  the  symmetric  values  given  in  equation  (1)  are 
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T(r.t)  -  0 
a(t)  •  a(t)  -  0 


(6) 


The  resulting  equation  and  boundary  conditions  cannot  be  satisfied  for 
arbitrary  a(t)  -  aQ  since  the  surface  tension  term  acts  as  an  inhouno- 
genelty  for  non-zero  aQ  and  finite  and  it  cannot  be  balanced  by  a 
toroidal  field  satisfying  all  three  constraints.  Thus,  to  conform  to 
(6)  the  degenerate  case  is  required 

T(r,t)  -  0  } 

f  t  <  0  (7) 

a(r,t)  •  0  > 

However,  without  some  additional  disturbance  only  the  trivial 
case  of  zero  nonsphericities  is  invoked.  The  necessary  inhomogeneity 
may  be  introduced  by  simply  specifying  «n(0)  and  an(0).  Physically, 
this  seems  troublesome  since  a  discontinuity  in  shape  and/or  nonspherical 
velocity  is  necessary.  In  the  viscous  case  this  objection  Is  not  rele¬ 
vant  since  both  mathematically  and  physically  the  manner  by  which  these 
initial  values  arise  doe  not  Influence  subsequent  behavior.  Thus ,  even 
if  a  discontinuity  in  aQ(t)  at  tQ  is  specified  which  requires  a  physi¬ 
cally  unreasonable  singularity  in  the  acceleration  the  problem 

is  well-posed  after  tQ.  These  same  arguments  cannot  be  forwarded  for 
a  viscoelastic  medium,  although  conditions  (6)  do  allow  values  to  be 
assigned  to  the  elastic  integrals.  For  the  viscoelastic  fluid  the 
manner  by  which  the  condition  at  time  t  •  0  arise  la  relevant  to 
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subsequent  (t  >  0)  behavior. 

Results  of  simulations  which  incorporate  the  pre-zero  static 
conditions  (7)  show  only  a  small  difference  between  a  viscous  fluid 
and  one  exhibiting  viscoelasticity  modelled  by  equation  (III. 9) 
with  X  >  10" 5  seconds.  (See  Results.  Section  VI.)  This  is  consist¬ 
ent  with  the  expectation  that  viscoelasticity,  modelled  in  this  man¬ 
ner,  will  be  cumulative  and  the  observation  that  in  the  present  case 
the  duration  of  the  nonspherical  bubble  dynamics  is  insufficient  to 
allow  stress  accumulation. 

The  difference  between  viscous  and  viscoelastic  responses  in 
this  analysis  results  from  a  difference  in  stresses  symbolized  by 

A (stresses)  • 

*  viscoelastic  stresses  -  viscous  stresses  (8) 

fT  e  (H’.t')  e  (H.x) 

-  f  (R.H)  J  N(t-t')  f(R.  >H»)  dT'  “  vf  (R,H)^R,hT 


where 

eg  -  strain  rate  in  viscoelastic  medium 
ey  •  strain  rate  in  viscous  fluid 

Thus,  the  difference  in  stresses  can  be  viewed  as  resulting  from  three 
highly  inter-related  and  Interdependent  components:  the  difference 
between  the  corresponding  instantaneous  strain  rates  e#  and  ey,  a  con¬ 
tribution  caused  by  the  0(c°)  convection  represented  by  f(R,H)  which 
might  be  called  "geometric''  and  finally,  the  existence  of  s  fundamental 
constitutive  distinction  between  the  "memory  functions”  N(t)  and  v  «(t). 
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Initially  neglecting  the  flret  two  contributions  so  that 


e  «  e,  ■  e  <  e 
e  v  — 


and 


fCK.H)  j  N(t-t-)  If 

— «D 

T 

S  J  N(t-t’)  ee(H',T')d  *  . 


(9) 


a  first  estimate  of  the  overall  stress  difference  can  be  derived 


T 

£  (stress)  -  |  (y)exp(yp— )e(T')dT' 
±  v  e«ax  I1-cxp(I)1 


(10) 


s  (y)  v  e _  for  small  (t/A) 


In  evaluating  the  integral  no  contribution  for  all  time  previous  to 
t  •  0  is  allowed  since  the  set  of  initial  conditions.  Equations  (6)  & 
(7),  specify  no  stress  field  for  that  period  of  time,  which  means 


T 

j  (£)exp(£2l)e(T')dT’ 

T 

-  j  (y)exp  (£~-)e  (t  *  )dr ' 


(11) 
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This  approximation  la  useful  tor  preliminary  at tamp t a  to  evaluate 
paramatar  sensitivity.  Ona  elreuastanea  In  which  thin  result  mould  ba 
Inaccurate  la  that  of  a  largo  temporal  variation  in  the  geometric  fac¬ 
tor  f(lil),  rendering  equation  (5)  Invalid.  However,  while  the  result- 
lng  atraaa  difference  may  be  large,  the  variation  In  f(E,H)  obtalne  for 
large  changes  In  1(t)  which  also  corresponds  to  an  lnereaae  In  Inertial 
effects  likely  to  decrease  the  net  effect  of  etreea  on  the  amplitudes 

s_(t).  Another  altuatlon  which  must  ha  Investigated  alnce  it  would 
n 

Increase  viscoelastic  effects  Is  any  aat  of  conditions  which  would  make 
Equation  (11)  Inaccurate.  The  time  period  before  t  ■  0  repreaents  the 
history  of  the  fluid  which  can  be  important  for  flulda  with  memory,  and 
a  self-consistent  scheme  to  describe  this  domain  should  be  found. 

A  first  attempt  to  Initiate  the  stress  integral  might  be  to  spec¬ 
ify  a  non-zero  constant  amplitude  ao  In  addition  to  the  conditions 
stated  in  (6).  However,  a  static  Inhomogeneous  solution  to  the  system 
does  not  exist.  Apparently  the  transient  problem  must  be  treated  to 
Introduce  nonsphericities  and  an  asymmetric  stress  field  prior  to  the 
Initiation  of  spherical  dynamics.  The  analysis  will  now  divide  the 
problem  solution  Into  two  time  intervals.  One  of  these  Intervals,  the 
second,  corresponds  to  the  full  problem;  after  time  t  •  t  ■  0,  the  system 
la  solved  for  aQ(T),  while  the  spherical  radius,  K(t),  la  allowed  to  vary 

In  time.  Before  tine  aero,  in  the  first  time  soma,  B(t)  •  1  ,  a  constant, 

o 

and  a  transient  scheme  to  Introduce  and  describe  asyMctrles  in  the  system 
must  be  devised.  This  proposed  theoretical  scheme  should  ba  mathematically 
tractable  and  also  allow  a  reasonable  physical  Interpretation.  Zf  possi¬ 
ble,  this  physical  significance  should  ba  directly  translatable  to  exper¬ 
imentally  achievable  conditions. 


Sine*  R(t)  »  Ro  and  i(t)  •  R(t)  -  0  for  all  tin*  t  ■  t  <  0  the 
complications  Introduced  by  a  La gran gl an  framework  ara  not  praaant  and 
their  (r,t)  coordinate  ayaten  can  be  choaen.  The  toroidal  field  equation 


then  bee owes 


t 

|  H(t-t»)[^-^f^Tn]dt» 


(12) 


If  a  separable  aolution  la  assumed  to  exist 


T„<r.t)  -  H  <1  X,,^  Pn>J(r»  03) 

where  6  . (t)  are  function*  only  of  tine,  p_  .  only  of  r  and  y  .  are 

n,j  n,j  “O 

constants,  then 


pn,j(r)  *  *n,J  (_df^  "  “ fT1  PntJl 

a*) 

d0  1  r* 

and 

Vl^“J*(t't,)en.J(t’)dt’ 

(15) 

where  R  .  la  a  particular  constant  of  separation.  A  steady  state  con- 
n»j 

tributlon  to  the  solution  la 


eo(t)  -  l 

ljs(r)  -  X^f-)®*1  «■  X8*-  <^)° 


(16) 


since  this  part  of  the  velocity  field  la  bounded  as  r*«  ^ 
equation  (111.30) 


0.  Proa 
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The  temporal  equation  haa  solutions  of  the  fora 


6j(t)  -  exp  l(~  ± 

J 


(18) 


which  means  that  the  separation  constant  is  complex 


V 


(f)/|T+ir+1,,31 

t  +  1“3 


(19) 


for  the  relaxation  function  N(t).  By  requiring  the  solution  to  be 
bounded  as  t  -*■  the  decay  time  is  restricted  to  nonnegative  values. 

The  spatial  equation  (14)  can  be  transformed  into  a  version  of 
Bessel's  equation.  For  n  *  2»  the  two  independent  modified  spherical 
Bessel  functions  relevant  to  the  solution  are 


V*>  •  -  h  *  ;>** 

(20) 

F2(x)  -  (j*r  +  +  £>•"* 


for  the  transformations 


x 


p(r)  -  r  F(r) 


Fj(x)  is  not  bounded  as  r  •,  so  only  Fj(x)  contributes  to  the  solu¬ 
tion.  The  spatial  integral  a^Ct)  consist  of  teree  of  the  form 
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°2.j(t)  ■  Vc)  *2.j  ('7> 


exp( 


,  dS 


(21) 


(22) 


when  D.  -  ✓SfT  /R 
J  j  o 

The  tine  dependent  portion  of  this  toroidal  function  is  an  exponen¬ 
tially  growing  oscillation  of  frequency  u>j  which  occurs  between  t  ■  -<*> 
and  t  ■  0.  The  time  constant  and  frequency  of  6^ (t)  can  he  specified 
Independent  of  any  variables  or  parameters  of  the  physical  system  as  cur¬ 
rently  described.  This  will  fix  the  separation  constant  Kj  which,  in  turn, 
determines  the  spatial  function  o^(r).  When  an  individual  transient  mode 
is  excited,  it  is  the  form  of  the  time  dependence  of  all  the  dynamics. 

Thus,  if  a  physical  mechanism  or  apparatus  could  be  visualized  to  drive 
some  event  with  form  0^ (t),  all  asymmetries  would  be  similarly  driven.  In 
0(e°)  dynamics.  Initiation  was  accomplished  by  changing  the  far-fleld  con¬ 
ditions  through  the  ambient  pressure.  An  initiation  of  0(e1)  dynamics 
through  far-field  conditions  would  be  ideal. 

The  oscillatory  time-dependence  is  reminiscent  of  rheological 
taatlng  which  employs  oscillatory  motions.  In  this  analysis,  far  field 
flow  conditions  are  governed  by  an(t)  in  *a(r,6,+,t)  (see  equations  III. 29 
and  III. 30).  The  value  of  a0(t)  in  equation  (III. 30)  results  when  the 
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velocity  field  vj  is  required  to  decay  as  r  increases.  A  way  to  alter 
that  requirement  is  to  add  a  time  dependent  term  to  the  old  aQ(t)  expres- 


(t)  -  an(t)  -  I  A j  e^t)  -  on(t)  -  l  Aj(t) 


which  is  formed  with  the  "appropriate"  time  dependence  6^(t).  The 
velocity  field  as  r  -*■  •  is 


yl  -  l  V  [Aj  e.j(t)rn  yJJ  (e,*)] 


For  the  particular  case  n  ■  2 


(yi)r  -*■  2Aj  e^(t)r  Y? 


(yi)e  -*■  ^  ®j<t)r  ||* 


Cyi)^  -*•  Aj  BjCt)  8me 


When  m  ■  0,  this  respresents  a  three-dimensional  elongational  flow  with 
extension  along  the  Z-  or  6  ■  0  -  axis,  and  time-dependent  elongation  rate 
^see  Figure  6a).  The  m  •  1  and  m  ■  2  cases,  with  n  ■  2,  are 
two-dimensional  elongational  flows  with  rates  Aj  9j (t)/6  In  the  y-s  and 
x-y  planes,  respectively.  (See  Figures  6b  and  6c.)  These  latter  two 
flows  are  particularly  convenient  since  they  can  very  nearly  be  achieved 
experimentally  in  the  Taylor  four-roller  apparatus. 


vtjUKWB 
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The  shapes  of  the  bubbles  created  by  these  flows  are  given  by 
(III. 12),  i.e.  the  deviation  fron  the  sphere  Is  given  by  (8,4)  with 
amplitude  a  (t).  Representative  shapes  are  shown  In  Figures  7(a)  and 
7(b),  8  and  9.  Figure  10  shows  the  streamlines  of  the  actual  external 
flow  around  a  spherical  body. 

The  new  far  field  condition  must  be  applied  to  the  boundary  condi- 

firm 

tlons.  All  a  (t)  must  be  replaced  by  a  (t).  The  flow  term  will  also 
n  n 

appear  in  the  pressure  term  Pv(r,t),  Equation  (III. 44)  since  it  arises 
from  the  equation  of  motion  which  employs  the  velocity,  not  vorticity 
and  the  velocity  field  Includes  the  influence  of  the  flow.  The  tangen¬ 
tial  stress  condition  for  time  t  <  0  requires  the  value  of  toroidal 
function  to  be 


Tj(R.t)  -  Xj  6j(t)RoDj(3Dj2+  3Dj  +  l)exp(-  ±  ) 

3 


(2$) 


and  with  a  new  intermediate  function  (t) 


“  Vt)110  Aj  +  3  XJ  •xp(“  5  )Dj(5V  +  5dj  +  1)f 

3 


(27) 


that 


a2(t)  -  a, 


*o  c  V° 

♦  g-  2,  -T-*1 - 

J  C  +  i“3 


(28) 


Again  the  constant  displacement  a-tt  must  be  taro  as  in  Equation  (7).  The 
normal  stress  balance  than  completely  determines  the  solution  for  all 
times  t  <  0.  The  final  result  is  expressed  through  e  cascade  of  new 
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Figure  8.  Representative  Bubble  Shape  : 

and  External  Flow  for  YgCe,*) 

Stream  surfaces  are  shown  along  with 
bubble  shape. 


Figure  10.  Alteration  of  Streamlines  by  Sphere 

for 

without  sphere  *•  pure  extension 
with  sphere 
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coop lex  variable* 


6  -  y  ,  for  fluid  relaxation  time  X 


y.  -  -  ±  iwj 
Go/p 


v  + 


3+y, 


D 

/  °\  /8  v  .  4o 

yc  ■  (T)yc  -  <r)JB +  f^r 

R 


yD  "  “P('^K  B  > 


yBT  "  <B  +  y53  Co/i>+v(S+y,) 


yS  ‘  yDI3yBT+  f)V(/y.1 
o 


(29) 


with  the  final  result 


40Ro  yt  A1 

-  I-  no  +  (-3-^)  -5-}  -i 
yc  ys 


(30) 


Evaluation  of  the  velocity  and  stress  fields,  as  well  as  the  bubble 
shape  is  now  possible. 

Notice  that  Xj  Is  complex  and  varies  linearly  with  the  complex 
strain  rate  A^.  Large  stresses  and  nonsphericities  can  be  generated  by 
specifying  a  large  (Aj  ( ,  but  are  subject  to  restraints  on  the  validity 
of  the  viscoelastic  model,  as  discussed  earlier, and  by  experimental 
practicability. 


( 
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The  behavior  o£  the  external  flow  after  8pherical  dynamics  are 
initiated  at  time  t  ■  0  is  a  new  consideration.  If  ia  not  large 
compared  to  the  positive  time  period  of  interest,  then  the  elonga¬ 
tions  1  rate  will  grow  to  large  values.  Therefore,  at  time  t  •  0,  a 
new  damping  factor  might  be  assigned  to  the  flow  and  since  the  full 
equations  will  be  solved  numerically,  without  the  separation  of  varia¬ 
bles  in  equation  (13),  this  is  mathematically  allowable. 

The  Introduction  of  the  external  flow  terms  allows  the  problem 
to  be  physically  and  mathematically  well-posed  at  time  t  ■  t  ■  0.  All 
quantities  now  have  physically  meaningful,  unambiguous  values.  Look¬ 
ing  toward  the  expected  behavior  for  subsequent  times  (t  >0),  the 
same  intuitive  arguments  presented  as  expressions  (8) ,  (9)  and  (10) 
suggest  the  kind  of  parameter  sensitivity  to  be  expected;  any  combina¬ 
tion  of  parameters  which  correspond  to  Bmall  instantaneous  strain 
rates,  but  allow  accumulation  over  long  times,  e.g  small  and  t3  »  X. 
would  show  the  largest  elastic  changes. 
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V.  EXPERIMENTAL  PROGRAM 

1.  Introduction 

In  conjunction  with  the  theoretical  analysis  of  nonspherical 
bubbles  in  viscous  and  viscoelastic  fluids,  an  experimental  program 
was  begun  to  test  the  model  and  attempt  to  Isolate  the  isost  Important 
factors  in  real  cavitation  inhibition.  The  use  of  dilute  polymer 
solutions  was  natural  in  view  of  the  initial  motivation  for  this  work, 
yet  added  complications  due  to  the  lack  of  rheological  measurements 
available  on  such  fluids.  Results  confirm  the  prediction  that  inertia 
would  dominate  spherical  bubble  dynamics  and  also  demonstrated  the 
viability  of  our  technique  in  creating  nonspherical  cavities. 

The  criteria  for  the  experimental  system  were  very  straight" 
forward;  it  must  have  the  capability  to 

1)  create  spherical  bubbles  which  remain  spherical 

2)  create  nonspherical  bubbles  reproducibly 

3)  allow  introduction  of  various  flows 

4)  record  these  bubble  sizes  and  shapes. 

Of  course  all  of  this  must  be  possible  in  various  fluids  with  some 
flexibility  with  regard  to  bubble  size. 

Previous  workers  had  used  numerous  techniques  to  study  single 
bubble  dynamics.  The  simplest  is  probably  the  Introduction  of  a 
cavity  by  inflation  using  a  device  such  as  a  syringe,  then  withdraw¬ 
ing  the  Instrument  snd  utilizing  surface  tension  to  seal  the  bubble. 

A  variation  on  the  same  method  depends  upon  buoyancy  to  pull  the  gas 
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volume  off  the  "Inflating  needle". 

The  technique  has  been  employed  with  success  by  Pearson  and 
Middleman  (1974,  1977,  1978)  to  generate  elongatlonal  flows  as  a  tool 
in  fluid  characterization.  However,  their  elongation  rates  of  order 
1  second**1  are  much  smaller  than  those  relevant  to  cavitation  bubble 
collapse,  where  the  time  scales  are  much  smaller  and  dissipative 
mechanisms  much  less  important  to  spherically  symmetric  motion.  As 
a  consequence  of  the  Reynolds  number  regime  into  which  Pearson  and 
Middleman's  work  can  be  classified,  useful  results  can  be  extracted 
without  regard  for  any  nonspherica-1  shape  uBlng  only  a  pressure  param¬ 
eter.  Thus,  the  presence  of  a  physical  apparatus,  which  Introduces 
asymmetries  into  the  bubble  environment  was  not  an  impediment  to  that 
work,  but  would  be  when  small  nonsphericities  are  important. 

One  way  to  remove  the  apparatus  from  the  imnedlate  vicinity  of 
the  bubble  has  been  to  allow  the  cavity  to  move  away  from  any  bubble 
forming  device  through  buoyancy  driven  motion.  Hassager  (1977,  1979), 
Coutenceau  and  Hajjam  (1981)  and  others  have  studied  both  viscous  and 
viscoelastic  effects  in  this  manner.  However  their  bubbles  are  usually 
large,  have  constant  volume,  and  the  first  order  effects  are  not  inertial. 
Once  again,  nonsphericities  arise,  but  it  is  the  steady  shape  of  the 
bubble  which  is  of  primary  interest. 

Experimental  events  with  dynamics  more  closely  resembling  cavi¬ 
tation  bubbles  can  be  generated  by  the  rapid  discharge  of  electrical 
potentials  of  the  kilovolt  range  across  electrodes  in  the  working 
fluid.  Gibson  and  Blake  as  well  as  Cbahlne  and  his  co-workers  have 
had  considerable  success  with  this  technique  In  studying  the  Inter- 
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action  of  the  bubble  with  various  controlled  asymmetrical  environments 
such  as  solid  walls  or  free  surfaces.  Their  configurations  confer  a 
planar  symmetry  on  their  experiments  which  appears  to  be  unimportant 
for  the  study  of  the  highly  deformed  shapes  which  arise  due  to  the 
presence  of  a  planar  boundary.  None  of  these  experiments  includes  the 
velocity  fields  which  are  present  in  flow  cavitation.  The  flow  pat¬ 
terns  might  well  be  altered  by  the  presence  of  velocities  in  the  regions 
of  high  shear  where  cavities  should  be  formed,  if  these  devices  were 
used  with  flow. 


2.  Apparatus 

A  technique  which  requires  no  apparatus  in  the  fluid  in  the 
region  where  bubble  forms,  has  been  pioneered  by  Lauterborn  (1981)  and 
his  co-workers.  A  laser  can  be  used  to  focus  electromagnetic  energy 
sufficiently  to  generate  cavities.  The  non-invaslve  character  of  this 
method  makes  it  the  most  promising  and  flexible  to  study  the  variety 
of  conditions  necessary  for  the  testing  of  the  model  system  analysis. 

The  apparatus  used  for  the  present  experiments  is  represented 
schematically  in  Figure  11  and  with  a  photo.  Figure  12.  The  major 
component,  a  ruby /glass  laser,  generates  maximum  output  of  1.5  J  in 
15-30  nsec.  The  effective  width  of  this  beam  as  detected  by  spot 
burning  tests  was  about  1  cm,  which  sised  the  focussing  lens  system. 
The  final  design  for  the  cavitation  chamber  without  flow  was  a 
parallelepiped  constructed  from  mMA  with  a  useful  Interior  which 
approximates  a  cube  11.5  cm  per  side.  The  wall  through  which  the 
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vtirure  11  Major  Components  of  Optical 
Figure  j  Cavitation  Apparatus 


Figure  12. 

Arrangement  of  Apparatus  Shoving: 

(A)  control  unit,  <B)  digital  time  delay, 
(C)  ruby  laser,  (D)  flash  lamp  and  diffuser 
(E)  test  cell,  and  (F)  Basselblad  camera 
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laser  bean  enters  the  tank  contained  a  piano-aspheric  lens,  mounted  flush 


with  the  interior  surface  and  properly  aised  to  focus  the  beam  to  a  crit¬ 
ical  volume. 

By  changing  the  single  lens  or  employing  a  series  of  two  lenses, 
the  focal  length  of  the  optics  could  be  altered.  This  was  useful  to 
evaluate  the  effect  of  the  proximity  of  the  lens  surface  which  acts  as  a 


solid  planar  wall  and  also  allowed  the  geometry  of  the  bean  focus  to  be 
changed.  This  becomes  particularly  important  aince  it  allows  generation 
of  nonspherlcal  bubbles  without  far-field  asyanetrles. 


The  bubble 


icorded  photographically  with  a  single  fr 


camera.  Since  the  time  scale  of  the  phenomena  is  tens  of  microseconds 
conventional  mechanical  shutters  are  not  adequate  for  the  resolution 
necessary.  Instead,  photo  exposure  was  accomplished  with  a  digital  time 
delay  of  rated  accuracy  ±0.1  ysec  triggering  a  flash  to  backlight  the 
image  which  is  transferred  through  an  open  camera  shutter.  High  speed 
instant  Polaroid  black  and  white  film  captured  images  magnified  four 
times  (4X)  by  the  lens  and  bellows  system. 

The  apparatus  listed  in  Table  4  allows  three  major  adjustments  in 
operating  conditions. 

1)  Laser  power  infringing  upon  lens  system.  Gross  adjustment  is 
made  by  attenuating  the  beam  using  the  schot  glees  filters.  Finer  ad¬ 
justments  are  made  by  edjusting  the  charging  voltage  across  the  capaci¬ 
tor  bank.  To  Insure  lasing  of  the  rod  a  charge  at  least  70S  of  the 
maximum  voltage  was  necessary.  This  translates  to  (.7)*  (1.5  J)  s  .755  J. 
since  energy  is  proportional  to  the  square  of  the  voltage. 

2)  Focal  length  L^.  and 

3)  Beam  convergence  angle  0g. 
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TABLE  4.  Equip— nt  Specifications 


Alignment  Laser 

Hughes  Model  4020 
120  —  Be-Ne 

Power  Laser 

Apollo  Laser,  Zac.  Modal  5H 
Ruby /Glass  Q-switch 
Power  1.5  J/S 
Pulse  Width  15-30  nsec 

Digital  Daisy  Generator 
BMC  Model  7010 

used  in  range  0-99999  t  1  psec 
Camera 

Hasselblad  Model  500C 

w /Zeiss  Planar  1:2.8,  f  -  80  am  lens  and  12  “extension 


Film 

Polaroid  667  high  speed  (ASA  600) 
black  and  white  Instant  development  film 

Flash 

EG&G  Electro  Optics  Model  85302 
with  Fresnel  diffuser 


Laser  Focussing  Lens 


focal  length 
(in  air) 

diameter 

Small 

12  — 

17  — 

Medium 

18 

24 

Large 

25 

32.5 

Diverging  Lens 

-40 

21 

Melles  Griot  f 


01  LAG  002 
01  LAG  005 
01  LAG  117 
01  LPK  013 


*Meutral  Density  Filters  (Schot  Glass) 

Thickness  Designation 


Approximate 

Absorbance 


Low  1  — 

Medium  2  — 

High  3  — 


EG5 

5% 

EG4 

10-20Z 

EGA 

SOX 

^Supplied  by  Mr.  Edwin  To Inass,  Princeton  University  PPL 
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These  are  two  geometric  parameters  defining  the  geometric  charac¬ 
teristics  of  the  focussed  beam.  As  shown  in  Figure  13,  the  "focal 
length"  Is  defined  as  the  distance  from  the  planar  surface  of  the  flush- 
mounted  lens  to  the  point  of  highest  energy  concentration  which  is  also 
the  location  of  cavity-inception.  As  a  first  approximation  this  is 
simply  the  focal  length  in  air  of  the  optical  system  divided  by  the 
index  of  refraction  of  the  medium  (1.33  for  water).  The  convergence  is 
defined  as  that  angle  at  which  the  beam  is  funneled.  For  a  simple 
single  law  system,  these  two  parameters  are  not  independent  since  the 
diameter  of  the  incident  beam  if  fixed  (nominally  considered  DhyMTO  “  1  cm). 
Linear  optics  gives  the  simple  relationship 


Lf  "  (Dbeam/2)  008  ®f 


(1) 


If  a  second,  diverging,  lens  is  placed  in  series  with  the  plano¬ 
convex  converging  lens  it  will  Increase  the  beam  diameter  impinging 
upon  the  primary  lens  so  that  although  equation  (1)  is  still  valid, 
Dbeam  can  varied  Independently.  Of  course,  Lf  is  now  a  simple  func¬ 
tion  of  the  focal  length  of  both  lenses,  approximated  by 


_J _ !_  + 1 

1.33  L,  fi  f2 


tt> 


Lf ,  as  stated  previously,  can  be  varied  to  assess  *Vall  affects". 
Adjustment  in  6f  should  affect  bubble  shape.  Intuitively  there  should  be 
some  critical  energy  density  which  incites  cavity  formation,  large 
values  of  6f  will  localise  this  density  sufficiently  to  create  spherical 
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bubbles.  Small  values  of  8^  should  create  a  more  extended  volume  with 
this  critical  concentration  and  thus  elongated  bubbles  (see  Figure  13(b)) 


3.  Types  of  Trials 

Four  major  types  of  experimental  trials  were  performed: 

IA)  Spherical  Bubbles  in  H2O 

These  trials  were  performed  to  evaluate  system  reproducibility  and 
determine  the  distances  over  which  wall  effects  are  present. 

IB)  Spherical  Bubbles  in  Dilute  Solutions 

Although  previous  experiments  using  this  optical  cavitation  tech¬ 
nique  have  been  performed  in  820,  similar  trials  in  dilute  aqueous  poly¬ 
mer  solution  had  not  been  reported.  The  first  goal  of  these  experiments 
was  to  guarantee  the  method’s  reliability  in  such  solutions  since  scat¬ 
tering  and  differing  inception  characteristics  Slight  be  expected.  Once 
the  technique  was  established  the  radial  dynamics  in  polymer  solutions 
could  be  compared  to  the  results  in  pure  water. 

2)  Initially  Spherical  Bubbles  near  Planar  Surfaces  in  A)  H20 
and  B)  polymer  solutions. 

The  apparatus  affords  an  excellent  opportunity  to  investigate  the 
effect  of  the  proximity  of  a  solid  wall  on  the  shape  of  a  growing  and 
collapsing  bubble.  Although  presently,  no  theoretical  analysis  has 
been  developed  to  treat  the  case  for  a  viscoelastic  fluid,  qualitative 
comparison  can  be  made  in  experiments  very  similar  to  thoae  conducted 
by  ChAhine  (1981)  and  others. 
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3)  Nonspherlcal  Bubbles  in  an  "Infinite  Medium" 

A)  H2O  and  B)  Polymer  Solution 

The  mechanism  by  which  nonapherlclties  arise  is  the  lengthening 
of  the  "critical  volume"  described  previously.  Sequences  of  bubble 
photos  allow  a  detailed  comparison  of  the  evolution  of  experimental 
shape  as  coupled  with  the  growth  and  collapse  of  the  cavity  and  theo¬ 
retical  predictions.  This  can  be  carried  out  both  for  the  pure  H20 
and  various  polymer  solutions. 

4)  Initially  Nonspherlcal  Bubbles  near  Planar  Surfaces 

One  possible  cause  for  differences  between  experimentally  produced 
bubbles  near  solid  walls  and  real  cavitation  bubbles  is  a  "competition" 
between  jet  formation  and  any  other  nonspherlcal  inodes  which  would 
reduce  the  energy  available  to  the  re-entering  Jets  (see  Discussion). 

A  first  attempt  can  be  made  to  evaluate  this  possibility  by  producing 
nonspherlcal  bubbles  near  solid  walls  with  different  orientations  for 
that  wall. 


4.  Fluids 

(R) 

Distilled  water  filtered  through  3.0  ym  Micropore  filters  was 
employed  in  most  trials.  The  effect  of  boiling  to  degas  the  water  was 
also  tested.  Preliminary  trials  showed  that  the  distilled,  filtered 
water  was  sufficiently  clean  to  eliminate  stray  nucleatlon  and  allow 
reproducible  trials. 

The  polymer  solutions  were  prepared  from  polyacrylamide  at 
M.V  •v  1.2  x  10*  obtained  from  the  American  Cyanamid  Corporation.  The 


polymer  samples  were  approximately  20Z  hydrolyzed.  Solutions  by  the  slow 
addition  of  a  weighed  amount  of  polymer-to  a  measured  volume  of  distilled* 
deionized  H2O,  then  tumbled  on  rollers  to  insure  homogeneity  for  at  least 
twelve  (12)  hours.  Initial  trials  were  performed  for  solutions  with 
concentrations  of  approximately  500  ppm  since  these  definitely  exhibit 
cavitation  inhibition  in  water  tunnel  tests. 


5.  Photo  Analysis 

No  sophisticated  procedures  were  available  for  data  analysis. 
Spherical  bubble  radius  was  measured  directly  from  the  photographs  with 
an  accuracy  better  than  ±  0.05  cm/4.  The  sphericity  was  evaluated  by 
overlaying  circles  and  checking  for  visible  deviation  from  that  shape. 

NonBpherlcal  measurement  was  leas  straightforward  and  multiple 
techniques  were  evaluated.  Each  involved  projection  of  the  photograph 
onto  a  circular  polar  coordinate  grid  and  the  assumption  of  axial  symmetry 
about  the  center  line  of  the  laser  radiation.  These  analyses  also  re¬ 
quire  an  assumption  of  the  surface  harmonics  Legendre  polynomials 
which  contribute  to  the  shape. 

Appendix  D  contains  sample  data  and  more  detail  on  procedures  chosen. 
The  final  results  shewed  dlfferehces  in  magnitude*  but  not  in  trends  be¬ 
tween  the  different  numerical  fits  attempted  so  the  simplest  —  estimating 
R  and  a^  by  measuring  the  major  and  minor  axes  of  an  image  —  was  employed. 


VI .  RESULTS 


1.  Theoretical  Spherical  Profiles 

In  order  to  generate  nonspherical  results,  appropriate  radius 
profiles  to  drive  the  asymmetries  are  necessary.  Previous  workers 
have  concluded  that  dilute  polymer  solution  rheology  will  not  cause 
experimentally  detectable  changes  in  the  overall  size  of  a  growing 
or  collapsing  bubble.  Stil,  an  0(e°)  model  similar  to  those  of 
Fogler  and  Goddard  (1970,  1971),  Ting  (1975),  Too  and  Ban  (1982)  and 
others,  which  incorporates  viscoelasticity,  but  employs  the  parameter 
values  of  the  present  study  (see  Table  3),  is  useful  to  compare  visco¬ 
elastic  change  in  this  geometry  to  changes  in  nonspherical  behavior. 

As  presented  in  Section  III,  equation  (4),  an  equation  similar 
to  Rayleigh's  original  result  can  be  used  to  generate  bubble  profiles. 
The  pressure  initiated  dynamics  are  completely  specified  by  the  dimen¬ 
sionless  parameters  of  Table  5 


TABLE  5.  Dimensionless  Parameters 
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where  R  «  -= — *  10“ 3  cm.  Is  the  initial  equilibrium  bubble  size 

O  r  — r 
V  a 

for  the  present  pressures  and  surface  tension.  Also,  the  Reynolds 

Number  as  defined  here  employs  a  ratio  of  the  length  and  time  scales 

to  define  the  velocity  scale  since  no  other  velocity  scale  emerges. 

The  time  scale,  t  ,  must  still  be  defined.  An  6(c°)  time  scale  might 
o 

Involve  P*  since  it  determines  the  dynamics,  however,  from  the  bubble 
growth,  collapse  and  constant  radius  cases,  no  single  meaningful  time 
scale  emerges.  A  time  scale  independent  of  the  particular  dynamics  is 
the  surface  tension  time  scale 


which  relates  to  the  period  of  an  undamped  oscillation  in  both  the 
0(c°)  spherical  or  OCc1)  asymmetric  modes. 

Results  for  R(t)  vs.  t  for  varying  ir  are  shown  in  Figures  14  and 
15.  For  positive  values  of  tr,  the  bubble  grows  rapidly  until  the 
velocity  R(t)  approaches  a  constant.  Then  is  no  change  for  any  value 
of  X  visible  on  the  scale  of  the  graph  and  numerical  values  show 
changes  of  less  than  0.1Z.  As  v  decreases  and  becomes  negative, 
collapse  is  described.  Once  again,  in  Figure  15,  their  is  no  visible 
change  due  to  variation  in  X,  but  f(t)  does  not  approach  a  constant. 
The  acceleration  R(t)  appears  to  continually  Increase  in  magnitude, 
this  is  due  to  the  surface  tension  term  in  the  R  equation. 


Model  Nonspherlcal 


cs  without  External  Flow 


The  initiation  scheme  of  Section  IV,  with  quiescent  ambient  condi¬ 
tions,  requirea  only  a  choice  of  a  particular  harmonic,  n,  and  an(0) 
and  a^tO)  beyond  the  parameters  of  spherical  flow.  The  problem  which 
allows  the  mathematically  consistent  use  of  an  initial  T(r,0)  function 


which  is  identically  sero  is  a  (0)  •  a  and  a  (0)  •  0  from  equation 

n  on 

(111.32). 

For  x  ■  0,  a  bubble  of  constant  radius,  R  ,  the  nonlinear  nature 

o 

of  the  governing  equations  is  removed  and  Laplace  transform  techniques 
have  been  used  by  other  workers  to  generate  a*(t)  vs  t  profiles.  The 
in tegro-d inferential  character  of  the  system  is  also  removed  for  a 
purely  viscous  fluid  and  physically,  the  viscosity  is  expected  to 
damp  any  oscillations.  This  expectation  is  confirmed  in  Figure  16 
where  a  30-fold  Increase  in  viscosity  is  shown  to  damp  the  asynmetrles 
greatly. 

Since  viscoelasticity  restores  the  time-integration  terms  to  the 
model,  the  possibility  of  somewhat  different  behavior  exists  even  for 
»  ■  0.  However,  as  discussed  previously,  the  absence  of  bubble  his¬ 
tory  to  initiate  these  time-integrals  at  some  non-sero  value,  in  con¬ 


junction  with  the  short  times  over  which  the  non-sphericities  persist 
(10** s  seconds ),aerves  to  allow  little  "elastic  accumulation"  and  re¬ 
sults  in  viscoelastic  behavior  which  differ  little  from  those  which 
would  occur  with  a  slight  change  in  viscosity.  This  is  shown  In  Figure 
17.  However,  the  fine  details  of  the  difference  cannot  be  reproduced  by 
a  simple  change  in  viscosity.  Khan  an  elastic  element  is  introduced, 
the  damping  of  the  resulting  amplitude  wavs  dsersasss,  which  would  occur 
if  viscosity  decreased,  hut  the  period  of  the  oscillation  lacrosses 


Figure  17.  Nonspherical  Amplitude  vs.  Time  -  Viscoelasticity 

No  external  flow 

aD  =  0  (see  Table  3) 


which  would  result  for  e  wore  viscous  fluid.  These  results  agree  with 
those  of  Inge  and  Bark  (1981)  who  generated  results  using  techniques 
suitable  only  to  R(t)  -  0. 
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Cases  of  wore  interest  involve  non-zero  dimensionless  pressure 
change*  ir,  values.  Stability  analyses  predict  that  bubble  growth  will 
restrict  the  viscous  nonsphericities  causing  the  bubble  to  remain  nearly 
spherical.  Typical  results  presented  in  Figure  18  show  this  to  be  true 
for  the  viscous  fluid  which  has  nonsphericities  which  become  overdamped. 
However,  the  viscoelastic  fluid  in  Figure  18  exhibits  somewhat  different 
behavior,  here,  as  R(t)  becomes  almost  constant*  the  shape  of  the  bub¬ 
ble  ,  as  measured  by  aa(t)/R(t)  also  appears  to  approach  a  constant. 

Earlier  stability  analyses  suggest  that  conditions  which  include 
bubble  collapse  may  result  in  nonsphericities  which  grow  without  bound. 
Figures  19  and  20  show  typical  profiles  which  confirm  this  prediction. 
Good  numerical  convergence  was  achieved  even  for  such  large  amplitude 
values.  If  the  linearized  system  remains  valid  such  bubbles  would 
exhibit  breakup  or  toroidal  ring  formation  and  this  mechanism  for  the 
dissipation  of  inertial  and  surface  energy  requires  consideration. 
Experimental  trials  are  necessary  to  determine  if  these  predictions  are 
even  qualitatively  correct. 

3.  Experiments  with  Ho  External  Flow 

The  initial  trials  with  the  experimental  apparatus  were  to  guaran¬ 
tee  the  ability  of  the  syatam  to  generate  bubbles*  first  in  water  and 
then  in  polymer  solution,  and  of  the  photographic  system  to  capture 
images  of  those  cavities.  The  first  successful  trials  resulted  In 
photos ,  typical  saw  las  of  Which  are  shown  in  Figure  21.  Mots  that 
these  bubbles  are  quite  large,  with  a  radius  greater  than  0.S0  cm  at 
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Figure  19.  Nonspherical  Amplitude  vs.  Time- 
Collapse  I 

No  external  flow 

a0  »  0 

Collapse  shown 
(see  Table  3) 


Figure  20.  Honspherical  Amplitude  we.  Tim# 
•  Collapse  II 

Ho  external  flow 

Collapse  profile  shown 
(see  Ikblt  9) 
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1*050  psec.  They  arc  created  with  canter  in  approximately  1.50, cm 
fro*  the  closest  asymmetry.  The  solid  "wall"  of  the  planar  face  of 
the  plexiglaa  lens  influences  the  bubble  In  all  dynamic  cycles  other 
than  initial  growth.  The  reproducibility  of  these  bubbles  is  excel¬ 
lent  with  a  superinposition  of  the  images  usually  possible  even  includ¬ 
ing  some  of  the  fine  structure  of  the  bubble.  There  are  occasional 
"gliches",  bubbles  which  do  not  reproduce  trials  for  identical  condi¬ 
tions.  these  appear  to  have  a  number  of  possible  causes: 

A)  Lack  of  laser  reproducibility.  A  qualitative  indicator  of 
the  laser  pulse  characteristics  can  be  obtained  by  looking  at  the 
bright  flash  at  the  original  bubble  center  which  is  the  overexposure  of 
the  photo  plate  caused  by  the  laser  luminescence.  Although  these  trials 
have  been  selected  for  their  reproducibility*  some  of  the  variation 
possible  is  indicated  by  a  comparison  of  the  spot  in  the  frame  at 

525  psec  and  that  at  1570  usee.  The  latter  event  shows  a  secondary 
brightness  separated  from  the  main.  Such  differences  usually  corre¬ 
lated  well  with  unexpected  differences  in  the  cavities. 

B)  Overly  Bapid  Bepetitlon  of  Trials.  Early  in  testing  it  be¬ 
came  apparent  that  reproducibility  suffered  when  an  experiment  was 
repeated  as  quickly  as  possible*  limited  only  by  the  charging  cycle  of 
the  laser  capacitors  which  requires  5-10  seconds.  The  lack  of  repro¬ 
ducibility  was  particularly  evident  with  conditions  and  at  times  which 
showed  nonspherlcal  bubbles.  At  such  times  waves  were  also  sometimes 
noted  on  the  free  surface  of  the  fluid  indicating  that  the  fluid 
motions  Induced  by  the  first  bubble  were  not  completely  d«ped.  It  was 
found  that  an  Interval  of  60  seconds  between  trials  was  more  than  suf- 
flcisnt  to  elimlnats  this  effect. 
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A  pusxling  feature  of  th«  photographic  sequences  la  tha  recording 
of  no  bubble  image  for  tine  delays  of  500  yaec  or  leas.  Initial  hypothe- 
aas  of  a  bubble  Induction  tine  proved  to  be  unfounded.  The  delay  la 
cauaed  by  a  lag  between  digital  delay  triggering  and  the  Inner  firing 
algnal  in  the  circuitry  of  the  laater  power  supply. 

A.  Spherical  Bubbles 

Spherical  sequences  were  generated  with  a  criterion  for  repro¬ 
ducibility  baaed  on  the  measured  radius  between  two  trials  varying  by 
51  or  less.  The  viability  of  the  technique  is  demonstrated  in  Figure 
22  which  shows  data  for  spherical  radius  vs.  tine  in  distilled  HjO  and 
476  ppm  polyacrylamide  solutions.  Conditions  used  were 

medium  lens  *♦>  ffi  *  1.5  cm. 
medium  neutral  density  filter 
801  laser  charge 

Bach  of  the  polyacrylamide  data  bare  represent  the  mean  of  at 
leaet  three  trials  bounded  by  the  resulting  standard  deviation,  Is- 
produclblllty  is  excellent*  especially  considering  the  high  growth  and 
collapse  values  produced.  The  slightly  higher  collapse  rate  for  water 
nay  result  from  difference  in  laser  interactions  with  the  two  fluids. 

The  order  of  magnitude  analysis  of  Table  2  suggests  that  a  small  part 
of  the  difference  may  be  rheological*  since  they  occur  in  the  region 
of  low  velocity  where  Inertia  is  smallest.  This  would  be  eliminated 
in  smaller  bubbles  «hare  surface  tension  becomes  svea  more  important  > 
however  the  else  here  is  the  minimal  for  the  single  shot  technique 
where  t  10  usee  emerges  as  s  limit  on  reproducibility. 


•.  2000  f.2500 

Figure  21 :  Bubble  Sequence  in  Water 
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B.  Spherical  Mblw  near  Solid  Walla 

Although  no  theory  has  bean  developed  here  to  include  the  Influence 
of  solid  walls  or  free  surfaces  on  the  evaluation  of  bubble  shape,  quali¬ 
tative  results  were  of  interest  to  confirm  the  findings  of  previous  workers 
and  also  to  allow  comparison  with  future  trials  which  would  include  ini¬ 
tially  nonapherical  bubbles  near  the  planar  surface.  The  present  photos 
confirm  earlier  findings  that  a  solid  wall  induces  bulk  motion  of  the 
bubble  toward  the  wall  and  subsequent  jet  formation  as  collapse  proceeds, 
no  visible  difference  appeard  in  the  jets  in  the  different  fluids  (see 
Figure  23). 

C.  Nonapherical  Bubbles 

The  scheme  by  which  nonapherical  bubbles  can  be  produced  in  an  other¬ 
wise  quiescent  fluid  was  described  earlier.  This  possibility  has  been 
mentioned  by  Lauterborn  and  Ebellng  (1977),  but  no  quantitative  results 
exist  or  any  assurances  of  reproducibility.  The  problem  of  laser  "spot" 
reproducibility  mentioned  earlier  with  regard  to  spherical  bubbles  is 
even  more  crucial  here  since  the  distribution  of  laser  energy  along  a  line 
segment  is  needed  to  create  a  disturbance  of  email  bubbles  which  grow  and 
merge  into  one  nonspherlcal  cavity. 

Initial  trials  were  encouraging  since  nonspherlcal  bubbles  were 
seen  to  be  feasible  but  they  suffered  from  an  extreme  degree  of  random¬ 
ness  in  the  particular  distribution  of  initial  nudeatic  sites  and 
their  strength.  The  line  of  cylinder  which  bounded  the  nuclaatlon 
region  was  well  defined  but  the  distribution  of  bubble  formation  within 
the  region  was  sporadic  which  caused  reproducibility  to  suffer.  More 
power  was  necessary  to  create  reproducible  bubbles  due  to  the  increased 


Comparison  of  Experimentally  Determined  Nonspherlcal  Amplitudes 
to  Theoretically  Generated  Curves  1 

;  -  distilled  eater  Theory  -  a0  ■  5.30  x  10 

Laser :  97%,  no  filter 
Lens  separation:  6  cm. 
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length  of  the  optical  path  In  the  fluid  (4-6  cm)  which  served  to  absorb 
and  scatter  some  light  energy  and  to  the  intentionally  lengthened 
region  in  which  sufficient  energy  concentration  was  desired. 

This  power  requirement  resulted  in  narrow  range  of  power  in 
which  nonspherlcal  bubbles  could  be  generated.  The  lower  bound  was 
determined  by  the  factor  just  cited  and  corresponded  to  the  low  den¬ 
sity  filter  and  85Z  charge.  The  energy  of  this  pulse  is  nominally 
1.0  J.  The  maximum  is  restricted  by  the  laser,  1.5  J.  This  range  is 
considerably  narrower  than  the  corresponding  one  for  spherical  bubbles 
which  could  be  produced  by  .50  J  and  less.  The  resulting  variability 
in  nonspherlcal  bubble  slse  was  small  with  maximum  equivalent  radius 
between  .3  and  .6  cm.,  but  those  which  could  be  generated  were  clear 
and  reproducible. 

Complete  sequences  were  generated  and  analyzed  for  distilled 
water  at  two  power  settings  and  two  polyacrylamide  solutions  near 
500  ppm.  The  bubble  shape  and  size  were  analyzed  (see  Appendix  D) 
and  the  results  are  shown  in  Figures  24-27.  As  discussed  in  Appendix 
D.l,  no  best  fitting  procedure  could  be  established,  all  gave  compar¬ 
able  results ,  so  that  the  simplest ,  measuring  major  and  minor  axes 
and  assuming  a  shape  conforming  to  the  second  harmonic  was  used. 

Thus,  the  error  bars  in  the  experimental  data  result  from  three 
sources : 

1)  reproducibility  of  bubbles 

2)  photo  image  measurement 

3)  uncertainty  in  fit. 

The  relative  importance  of  these  factors  varied  from  image  to  image 


Figure  23:  Initially  Spherical 
Bubble  Collapsing  near  a 
Solid  Wall  -  Impinging 
Jet  Formation 
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Figure  26.  Comparison  of  Ixperlaeatally  Determined  Non 
spherical  Amplitudes  to  Theoretically  Generated 
Values  III 

Experiment  -  Laser:  17.81,  low  filter 
Leas  separatioa:  4.5  ca. 

Theory  -  a«  •  1.27  x  10”'  cm. 

&o  -  —4.0  x  10  ea./sec. 
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Figure  37.  Coapulwn  of  Experimentally  Determined  Nonsphevlcal 
Amplitudes  to  Theoretically  Generated  Values  IV 

Experiment  -  Laser:  07%,  no  filter  Theory  -  -  5.0  x  10'i  cm. 

Lens  separation:  6  cm.  40  -  -3.0  x  10*  cai./sec 
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within  •  mne>i  although  factor  2)  waa  usually  laaat  important.  Da* 
spite  thaaa  uncart aintiaa  a  aaooth  prograaaion  of  ahapa  waa  tha  raault. 

Tha  K(t)  n.  t  profile  waa  awoothad  by  a  boat  fit  to  a  tourier 

aarlaa  whan  Initial  tine,  period  and  tha  contribution  and  elgenvdua 

of  up  to  five  eigenfunctions  ware  varied.  Thaaa  fita  are  ahovn  aa 

tha  "center-line"  dotted  curve  in  Figures  24-27.  This  fit  waa  than 

differentiated  twice  for  uae  aa  the  forcing  function  in  the  aolutlon  of 

the  nonspherlcal  naplltude  equations.  The  changes  in  sign  of  the  radial 

acceleration  R(t)  in  the  fit  profile,  corresponding  to  inflections  in  the 

curves,  occur  because  of  the  fitting  procedure  and  the  varying  nuXber  of 

data  points  taken  at  different  tines.  Despite  the  resulting  physically 

unrealistic  values,  the  nonspherlcal  amplitude  profiles  driven  by  the  fit 

are  well-behaved.  The  curves  are  snooth  because  the  inflections  occur 

while  the  bubble  is  growing,  a  period  during  which  nonsphericities  ere 

stable  and  relatively  insensitive  to  the  details  of  the  radius  profile. 

The  initiation  of  the  nodal  required  an  initial  tine  tQ  at  which 

the  anplltude  a(t  )  had  been  neasured.  This  tine,  t  ,  was  chosen, 
o  o 

before  any  date-specific  nodelling  was  done,  by  taking  it  as  the  tine 

of  the  first  data  point.  The  first  data  point  corresponds  to  the 

first  bubble  photo  in  a  tine  sequence  which  appeara  as  one,  relatively 

aaooth.  cavity  and  usually  occurred  by  55 0  usee  or  50  usee  after  any 

bubble  is  visible.  Since  there  was  no  external  flow  and  notion  was 

present  for  less  than  50  usee  previous  to  t0*  bubble-induced  fluid 

history  was  Ignored;  the  stress  integrsls  begin  with  sero  value. 

Variation  of  *(t  )  and  fluid  paraaaters  was  than  perforasd  to  optldse 
o 

the  agreaaent  with  experlneatd  data. 

tatting  a(t0)  to  various  values  and  varying  fluid  paraaaters 
showed  that  viscoelastic  affects  warn  snail  as  would  be  expected 
fron  previous  nodal  tuna  for  situations  without  extend  flow  and  with 
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zero  itrui  integrals.  Thus  than  ara  two  important  da  grass  of  freedom, 
tbs  Initial  nlocity  a(tQ)  sad  tbs  viscosity  p.  lbs  sbsanca  of  viaco- 
alsstlc  offsets  also  rsdoesa  tbs  sensitivity  of  tbs  agr slant  to  tbs 
initial  tins  tQ.  Changing  tQ  will  not  altar  tba  nodal  predictions  for 
a  given  u  and  corresponding  initial  conditions ,  it  can  only  change  tba 
number  of  data  points  to  be  fit  and  bars  that  number  is  already  maxi¬ 
mized. 

lbs  viscosity  (or  Reynolds  Dunbar)  values  which  gave  tba  good 
agreement  shown  by  the  solid  lines  vara  lever  (or  higher)  than  expected 
for  water.  They  were  between  0.1  and  0.3  cp.  Setting  the  viscosity  to 
1.0  cp  without  altering  tbs  initial  velocity  resulted  in  predictions  of 
overdamped  oscillations  during  bubble  growth;  the  amplitude  S2(t)  did 
not  cross  the  S2  ■  0  axis  while  radial  velocities  were  positive.  In¬ 
creasing  the  magnitude  of  the  initial  velocity  a(tQ)  and  the  viscosity 
did  serve  to  drive  the  amplitude  past  spherical  equilibrium,  but  the 
shape  was  altered,  sea  Figure  24,  and  fit  was  lass  good. 

The  most  encouraging  aspect  of  the  model  results  may  be  the  accu¬ 
rate  predictions  of  trends  even  whan  deformations  become  large  so  that 
the  linearised  analysis  should  break  down.  In  both  water  and  solution, 
the  final  frames  of  a  sequence,  as  tbs  bubble  volume  became  very  small, 
could  show  one  of  two  distinctly  different  "asymptotic”  behaviors.  One 
type  is  suggested  by  tbs  upper  left-hand  illustration  in  Figure  9;  the 
nonspherical  amplitude  is  positive  and  becomes  large  as  R  decreases  so 
that  the  cavity  is  extended  in  a  line  along  tbs  axis  of  symmetry.  The 
other  extreme  occurs  for  negative  amplitudes  and  gives  tbs  forms  shown 
in  the  lower  right  of  Figure  9*  In  experimental  photos  this  appears 
as  an  image  like  Figure  28(a)  and  becoams  a  vertical  line  as  R  ♦  0. 

Such  shapes  are  too  distorted  from  tbs  sphere  to  be  represented  by  an 
equivalent  radius  R  and  amplitude  at.,  and  predicted  by  theory* 
Still,  the  model  successfully  predicts  these  final  trends  which  is 
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extremely  useful  for  the  study  of  bubble  breakup. 

D.  Nonspherlcal  Bubbles  Near  Solid  Wells 

The  preaise  upon  which  this  work  Is  based  is  that  the  role  of 
velocity  fields  in  cavitation  dynamics  extends  beyond  the  creation  of 
the  dynamic  pressure  difference  necessary  for  nucleation,  that  it  alters 
the  shape  and  the  size,  i.e.  the  total  behavior,  of  the  subsequent  cav¬ 
ity.  Another  factor  of  practical  importance  in  this  behavior  has  been 
shown,  both  here  and  elsewhere,  to  be  the  proximity  of  a  solid  wall, 
but  this  effect  has  only  been  investigated  for  initially  spherical  bub¬ 
bles.  Trials  were  conducted  to  Investigate  any  Interplay  between  non- 
spherical  bubbles  and  solid  walls.  They  are  complicated  considerably 
by  the  introduction  of  a  direction  vector  for  the  bubble,  which  can 
simply  be  chosen  as  along  the  initial  axis  of  symmetry  for  these  bub¬ 
bles,  and  it  can  vary  in  orientation  with  respect  to  the  solid  wall 
which  can  be  represented  by  a  normal  vector.  Quantitative  comparison 
to  results  for  spherical  bubbles  without  detailed  theory  is  ambiguous. 
The  distance  of  the  initial  center  of  the  bubble  from  the  solid  wall 
scaled  on  a  maximum  bubble  radius  gives  a  good  variable  with  which  to 
correlate  results  (Chahlne,  1981)  for  a  spherical  bubble.  However  for 
asymmetrical  bubbles,  the  amplitude  of  any  Initial  nonsphericity  along 
with  the  orientation  of  the  direction  vectors  will  Introduce  much  more 
variability  into  the  system  even  for  identical  fluids. 

As  a  result  of  these  complications,  only  qualitative  results 
were  sought  in  these  exploratory  experiments.  Comparisons  between 
BaO  and  polymer  solutions  were  also  mads. 


Photographic  evidence  from  these  trials  lead  to  two  tentative 
conclusions : 

1)  A  strong  competition  does  exist  between  the  jet-forming  effect 
of  a  solid  wall  and  nonsphericities  which  are  separate  from  that  effect 
and  it  varies  with  orientation. 

2)  Ho  large  differences  between  fluids  occur. 

This  first  conclusion  is  supported  by  the  photo  sequences  of  Fig¬ 
ures  30  and  31  which  show  no  images  of  a  well-formed  jet  like  the  one 
in  Figure  23.  Figure  23  differs  from  Figure  30  and  31  in  that  the 
Initial  optical  cavity  in  the  former  was  spherical  while  Intentionally 
nonspherlcal  bubbles  were  initiated  in  the  latter.  Figures  30  and  31 
are  representative  of  images  taken  over  a  range  of  bubble  sizes  and 
wall  proximities.  All  showed  bubble  migration  toward  the  wall  upon  col¬ 
lapse,  but  none  showed  an  impinging  jet.  Figure  30  shows  the  results 
when  B  and  N  are  perpendicular  (see  Figure  28(a)). 

When  B  and  H  form  a  45°  angle  more  experimental  problems  develop 
since  a  bubble  which  is  too  highly  nonspherlcal  Initially,  and  close  to 
the  wall  will  contact  it.  Thus  the  location  of  bubble  center  is  more 
limited  than  for  the  previous  cases.  Here,  again  the  solid  wall  causes 
bubble  migration  toward  it  with  distortion  but  no  jet  is  apparent  (see 
Figure  31). 

Trials  with  other  orientation  angles,  wall-bubble  distances  and 
500  ppm  polyacrylamide  solutions  yielded  similar  results.  Of  course, 
such  trials  are  not  comprehensive,  especially  since  the  nonspherlcal 
bubbles  are  vary  limited  in  their  slse  and  shape  prof  ilea.  However 
the  contrast  in  the  case  of  jet  formation  for  spherical  bubbles  and  the 
absence  thereof  for  asynnetrlcal  bodies  is  striking. 


Figure  29:  Bubble -Wall  Orientation 


a.  B  1  N  b.  B  45°  N 


Figure  30:  Sequence  -BIN  I  psec. ,  4x) 


rw 


4.  Model  with  External  Flow 


Recall  the  aodel  which  generates  a  completely  self-consistent  fluid 
history  by  postulating  an  external  flow.  The  details  are  developed  in 
Section  lv  and  summarised  here  and  in  Figure  32.  In  the  seal-infinite 
tiae  domain  with  upper  bound  time  t  ■  t  -  0,  the  equations  governing 
asymmetric  dynamics  can  be  driven  by  a  mathematical  forcing  function 
which  corresponds  to  an  externally  imposed  flow  in  the  physical  realm. 
Mathematically,  any  such  forcing  function  can  be  linearly  decomposed  into 
terms  with  time  dependencies  of  the  form 


©r  j(t)  -  6j (t)  -  exp[(~  t  irn  )t], 

J 


t  <  0 


(XV. 18) 


Since  this  forcing  function  is  bounded  throughout  the  Interval  [-*,  0], 
the  time  constant  is  positive  (and  real).  When  specific  values  of 

and  the  frequency  are  chosen,  they  completely  determine  the  forced 

response  of  the  fluid  for  a  given  angular  mode  n  and  fluid  model.  The 
stresses  in  the  fluid  and  shape  of  the  bubble  are  also  determined.  The 
terms  remaining  in  the  resulting  flow  field  as  r  •  are 


y?  -  Re  l  l  A  6  .(t)7[rV{6,*)] 

n  j  n’’  n*J  n 

-  **  I  I  An  .(t)7lrV(e,t)) 

n  J  J 


(IV. 24) 


where  -  Aq  is  a  complex  constant  and  Aq  ^  (t)  •  A^(t)  is  a  complex 
function  of  time.  For  n  •  2,  as  shown  in  Equation  (IV. 25)  and  Figures 
6,  8  and  10,  this  spatial  variation  corresponds  to  two  and  three- 
dimensional  extensions!  flows  with  tlms-dependent  alongatlon  rates 
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Figure  32.  Conparlson  of  Flow  and  Ho  Flow  Condition. 


Expressions  shown  for  n  »  2,  hut  corresponding 
conditions  exist  for  ell  n  >  2. 


directly  proportional  to  A  (t). 

n 

At  tine  zero  the  external  flow  haa  generated  nonsphericities  in 
the  babble  of  radius  Rq  and  the  added  complication  of  a  time  dependent 
equivalent  radius,  R(t),  can  be  Introduced.  (See  Figure  32.)  The 
particular  solution  obtained  by  the  separation  technique  of  Section  IV 
no  longer  applies  when  R  chsnges  so  that  the  full,  numerical,  procedure 
of  Appendix  C  eust  be  utilised.  Within  this  procedure  the  external 
flow  is  no  longer  constrained  to  the  time  dependence  of  expression  (IV. 16). 
To  investigate  the  Influence  of  a  further  alteration  of  the  flow,  the 
amplitude  of  the  oscillation  is  allowed  to  decrease  with  time,  but  the 
frequency  u  was  retained.  This  amplitude  variation  is  characterized  by 
a  time  constant,  ^decay*  Various  mathematical  forms  of  the  decay  function 
were  tried,  e.g.  a  ramp  or  a  sigmoid,  but  the  results  were  insensitive  to 
the  particular  choice  so  the  most  common  form  of  decay  was  chosen,  an 
exponential,  s.t. 


A  (t)  -  A  exp(-  •2- - )  sin(w.T) 

n  “  'decay  3 

for  t  >  0 


Physically,  this  damping  was  Introduced  to  simulate  a  change  in  flow 
environment  for  the  cavity  and  surrounding  fluid,  e.g.  migration  of  the 
bubble  from  a  region  of  high  extension  to  one  experiencing  lower  rates. 

The  time  dependence  of  the  external  flow  is  characterised  by  four 
parameters:  the  complex  amplitude  constant  Aq,  the  frequency  u  and  two 
time  scales  £  and  Tjacay  which  apply  in  the  Intervals  [-*,0]  and  [0,«], 
respectively.  Dimensionless  groups  which  Incorporate  these  values  ere 
shown  in  Table  6.  The  scaling  of  the  amplitude  a^  on  the  magnitude 
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{Aji  results  from  equations  (IV. 28)  end  (IV.  30)  which  show  that  these 
two  quantities  are  directly  proportional  at  time  t  •  t  ■  0. 


TABLE  6.  Dimensionless  Groups  for  External  Flow 

*2 

T^TV*o  <n  "  2) 


flow  UT 


Tflow  transient  "  Tf.t, 


-  — £■ 


decay 


decay 


Sensitivity  of  the  final  results  to  Tflw  tran8ient  was  easily 
anticipated.  As  the  value  of  t.  decreased,  viscoelastic  changes  were 
also  reduced.  This  is  consistent  with  the  importance  assigned  to  fluid 
history  in  distinguishing  viscous  and  viscoelastic  dynamics.  Since  this 
work  is  concerned  with  viscoelastic  changes,  this  effect  was  minimised 
by  allowing 


lf.t. 


(8) 


Sensitivities  to  and  are  discussed  in  conjunction  with  the 

upcoming  results. 

The  simplest  flow  case  obtains  for  w  ■  0,  a  bubble  of  constant 
volume,  just  as  it  does  for  the  quiescent  fluid  analysis,  kesults  for 


103 


tvo  different  flow  frequencies  or  values  ere  shown  in  Figure  33. 

The  only  difference  between  (s)  end  (b)  is  Figure  33(e)  shows 

thet  xflow  ■  2.5  *  10s  represents  low  frequency  or  slow  flow  since  saell 
lnstentsneoue  strain  retes  whose  effects  cen  be  multiplied  by  visco¬ 
elastic  stress  integration  result.  Thus,  e  rather  large  difference  is 
seen  between  the  viscous  and  viscoelastic  fluid  responses  and  it  increases 
with  increasing  Deborah  number.  However,  as  decreases,  the  effect 

of  viscoelastic  properties  decreases  (see  Figure  33(b)).  The  curves  do 
not  oscillate  about  the  a  ■  0  axis  because  for  these  flow  frequencies 
and  phase,  the  flow  is  unidirectional  over  this  time  span,  constantly 
"pulling"  the  bubble  in  one  direction. 

Conditions  which  combine  flow  and  growth  show  no  qualitatively  new 
flow  effects.  The  growth  damps  the  oscillations  and  does  not  allow  any 
initial  viscoelastic  differences  to  be  magnified. 

In  Figures  34  and  35  results  are  presented  for  different  collapse 
profiles  with  T^ecay  ■  0.25.  The  damping  of  the  external  flow  served 
to  magnify  differences  between  fluid,  which  were  present  at  t  -  0. 

This  is  entirely  consistent  with  the  notion  that  the  viscoelastic  fluid 
will  "remember"  the  past  occurrence  of  the  external  flow  while  the  vis¬ 
cous  fluid  responds  only  to  Its  Instantaneous  presence. 


FLOW 


n  *o 


a 


Figure  33.  Effect  of  External  Flow  Frequency 

(see  Table  3) 


VII.  CONCLUSIONS  AND  DISCUSSION 


The  results  of  the  previous  section  were  presented  individually, 
hut  it  is  only  When  they  ere  taken  collectively  that  more  exciting  con¬ 
clusions  can  be  drawn  about  real  flow  cavitation  and  their  contribution 
to  the  understanding  of  cavitation  inhibition  can  be  assessed.  In 
Table  1,  the  scope  of  the  present  work  was  resolved  into  three  cate¬ 
gories.  The  results  obtained  here  suggest  that: 

1)  Without  any  external  flow,  fluid  rheology,  and  slight  visco¬ 
elasticity  in  particular  do  not  strongly  Influence  cavitation  bubble 
dynamics  for  A)  spherical  or  B)  nonspherical  bubbles.  The  latter  cavi¬ 
ties  do  display  non-negligible  effects. 

2)  An  external  flow  which  is  the  sole  source  of  system  asym¬ 
metries  will  influence  bubble  dynamics  and  do  so  through  fluid  rheology. 

3)  The  characteristics  of  bubble  collapse  near  a  solid  wall  are 
altered  by  nonsphericities  which  are  present  due  to  influences  other  than 
the  "solid  wall  effect",  e.g.  external  flow. 

Of  course,  these  conclusions  all  hold  only  for  the  particular 
theoretical  and  experimental  systems  tested.  The  fluid  have  viscosities 
on  the  order  of  1.  cp.  and  the  elastic  number.  Ei,  has  a  value  near 
unity.  These  values  may  underestimate  the  Influence  of  viscoelasticity 
on  real  dilute  solution  since  oscillatory  viscometric  measurements  for 
similar  solutions  give  zero  shear  viscosities  approaching  1  poise 
(Chang,  1975)  but  the  values  employed  do  model  the  experimental  results 
very  well  which  may  be  due  to  the  shear-thinning  behavior  of  the  fluid. 
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There  Is  a  fundamental  rheological  difference  between  the  non- 
ephericitiea  in  conclusion  1)  and  cone lua ion  2).  In  the  foraer  case, 
rheology  is  Involved  only  in  the  evolution  of  the  as yame tries,  not 
their  genesis  or  initiation.  The  fluid  properties  actually  serve  to 
generate  the  shape  and  stress  field  in  the  latter  case,  through  the 
external  flow.  Since  rheology  is  important  in  both  generation  and 
evolution,  (pre-re ro  and  post-zero  dynamics),  it  is  more  Important  for 
this  case. 

Combined  with  conclusion  2),  3)  can  be  seen  as  an  example  of  a 
situation  where  viscoelasticity  can  be  very  Important  in  cavitation 
near  solid  walls.  Any  flow  past  the  wall  will  serve  to  create  a  stress 
field  surrounding  the  cavity  which  will  vary  according  to  fluid  rheology. 
This  will  create  a  nonspherical  bubble,  and  these  flow  Induced  asym¬ 
metries  will  effect  the  jet  Induced  by  the  solid  wall  upon  bubble 
collapse. 

This  work  also  makes  it  clear  why  many  previous  model  systems, 
which  were  constructed  in  an  attempt  to  display  large  viscoelastic 
effects  comparable  to  the  cavitation  Inhibition  of  Figure  1,  failed. 

These  were  designed  to  simulate  cavitation  in  a  quiescent  fluid,  and 
thus  Ignore  the  importance  of  stress  history  and  flow. 


VIII.  THE  NEXT  STEPS 


1.  Theory 

Experiments  have  shown  the  expansion  technique  employed  here  to 
be  valid  for  surprisingly  large  departures  from  spherical  shape.  How¬ 
ever  this  first  order  procedure  is  probably  as  far  as  a  linearised 
model  can  be  extended.  The  next  step  along  similar  lines  in  theory 
would  incorporate  a  fully  nonlinear  constitutive  model.  Some  of  the 
relative  simplicity  of  the  present  work  might  appear  retainable  through 
the  introduction  of  a  second  expansion  parameter,  this  one  characteriz¬ 
ing  the  flow  field  and  not  the  bubble  shape.  Then,  a  second  order  expan¬ 
sion  in  this  new  parameter  would  be  formulated  so  that  nonlinear  rheo¬ 
logical  effects,  especially  normal  stressea  are  generated.  Unfortunately, 
there  are  at  least  two  definite  problems  with  this  approach.  The  second 
order  terms  will  eliminate  the  use  of  the  linear  independence  of  the 
spherical  harmonics.  The  angular  dependencies  will  no  longer  be  separate 
from  the  (r-t)  formulation  (see  Appendix  B).  Even  if  this  were  resolved, 
the  size  of  the  bubble,  governed  by  the  zeroth  order  equations,  is  not 
coupled  with  the  asymmetries  and  there  is  an  interaction  for  the  large 
deformationa  that  the  second  order  expansion  would  attempt  to  predict. 

The  next  useful  step  in  any  modelling  effort  is  treatment  of  the 
full  problem.  This  is  an  extremely  difficult  undertaking  since  it  is 
highly  nonlinear,  transient,  contains  an  undetermined  free  surface  and 
requlras  an  infinite  domain.  Since  the  development  would  aim  toward 
incorporating  a  general  flow  field  and  a  solid  wall,  no  symmetry  sim¬ 
plifications  emerge  and  the  full  three-dimensional  problem  appears 
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necessary.  Marker  and  cell  techniques  which  have  been  used  to  model 
bubbles  near  solid  vslls  without  extemsl  flow  night  be  applicable  for 
the  Newtonian  liquid  or  even  a  shear  thinning ,  but  Inelastic ,  fluid, 
f  However,  Incorporation  of  non-Newtonian,  elastic,  effects  probably  re¬ 

quires  a  finite  element  treatment,  techniques  for  which  are  still 
under  development. 

1 

2.  Experiments 

Many  of  the  experiments  proposed  here  are  in  various  stages  of 
development  at  Princeton  University. 

The  early  successes  of  techniques  to  generate  nonspherlcal  bub¬ 
bles  with  or  without  solid  walls  present  suggests  that  similar  experi¬ 
ments  be  attempted  on  other  liquids  and  on  the  same  fluids  described 
here,  but  under  different  conditions.  The  aim  of  such  changes  in  ex¬ 
perimental  conditions  would  be  to  broaden  the  range  of  bubble  sixes 
and  profiles  for  which  reproducible  trials  can  be  conducted.  Alternative 
fluids  include  a  glycerine /H2O  solvent  system  which  allows  variation  in 
viscosity  through  the  composition  ratio,  and  could  be  investigated  with 
and  without  polyacrylamide  solutes.  More  flexibility  in  useful  trials 
for  the  aqueous  system  might  be  achieved  by  locally  or  non-local ly  dying 
or  tinting  the  fluid  to  change  energy  absorption  by  the  fluid.  Such  a 
procedure  would  need  careful  testing  for  such  complications  as  asymmetric 
thermal  effects  which  would  alter  reproducibility. 

All  trials  should  be  viewed  with  particular  regard  as  to  how  and 
when  jets  form.  These  variations  might  be  measured  using  a  pressure 
transducer  on  the  solid  wall  to  record  the  impact  of  «**  jet  which  forms. 


More  sophisticated  photographic  techniques  nay  also  allow  measurement  of 
jet  velocities. 

Flow  experiments  also  need  to  be  performed.  The  best  flows  would 
be  the  two-dimensional  extension  approximated  by  a  four-mill  or  four- 
roller  apparatus  or  two  Impinging  sheet-shaped  jets.  However  any  well- 
characterized  flow  which  imparts  stress  history  to  the  fluid  would  give 


useful  results. 
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APPENDIX  A:  o(e)  Equation* 


A.l.  Toroidal  Field  Equation 

The  full  velocity  field,  to  0(e),  is  given  by 

V  ■  v  +  ev  +  ev 
~  ~  o  p  ~v 


where  the  (e°)  field,  in  apherical  coordinate*,  i* 


<Ve  “  <*•>♦  "  0 


RR2 


<Yo>r  "  ? 


The  potential  contribution  derived  by  Plesset  is 

v  •  y  v  ♦ 

-p  £  -  p,n 

,  _n+2  • 

4  -  -  ~tt  ”T7T  la  +  2a  Y^(e,*) 

’p.n  n+1  ^n+1  n  n  R  n 


Rheological  contributions  are  contained  in  vy  term  which 
Equation  (111.28)  as 

3B 


<Vr  - 1 -  >rK 


B  3T* 

<V6  - 1  <-  r> r 
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(111.16) 


(III.l) 


(111.19) 


obtains  from 


(A.1) 


A-2 

where  B  (r.t)  ie  defined  as  the  function  containing  the  redial  and  tem- 
n 

poral  dependence  of  the  Integral  defined  by  equation  (III >29)  with 
the  external  flow  contribution  of  equation 

♦  (r,0,$,t)  -  B  (r.t)Y®(6,*)  (A-2) 

11  U  u 

Bn(r,t)  ■  [oft(t)  -  An  (t) 

+  jSrf  s‘n  V-t)dslr” 

R 

+  <Si  ^  !«■<*>  *  A„ 

(A.  3) 

+  &r  f^1  *'*H> 

R 

-  X <r,t)rn  ♦  Z  (r,t)r"(n+1) 
n  n 

This  equation  also  serves  to  define  X^  and  Z^,  with  °n(t)  given  by 
equation  (III . 30) . 

The  goal  of  this  first  derivation  is  an  analogue  to  the  vorticity 
equation 

lE+V'Vtf-w'Vy  -  Vx (stress)  (III. 22) 

9k  “*  m 

since  the  simultaneous  evaluation  of  the  pressure  field  and  toroidal 
function  T  is  not  necessary  nor  is  the  potential  flow  field.  This 


wit  be  accomplished  with  special  attention  to  the  Lagrangian  and  non- 
instantaneous  nature  of  the  stress  expression.  The  left-hand  side  (LHS) 
of  the  equation,  containing  the  terms  which  represent  inertial  effects 
is  not  altered  by  the  particular  constitutive  relation  employed,  hut  does 
need  to  be  expressed  in  the  proper  coordinate  system. 

The  components  of  the  vortlcity  vector  in  Eulerlan  coordinates  are 
first  order  in  e. 


(111.21) 
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(A.  4) 


The  (LHS)r,  the  radial-component  equation  of  the  left-hand  side  of 
the  vortlcity  equation,  has  terms  which  are  all  Identically  sero.  At 
first  order  in  e,  the  0  and  4  equations  vary  only  in  their  angular  de¬ 
pendencies 


<LHS)0  (LHS)^ 

1  3Y  3Y 

sin©  34  30 
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.  I  {  — —  +  1_  [a(£)2  v  1} 
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This  term  can  be  expressed  in  Lagrangian  (B,t)  coordinates,  defined  by 


equation  (111.38)  as 


3T  (H,t) 
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-  2  R*R 


T  (H,t) 


L» 


-} 


(A.  6) 


where 


L  -  (R*  +  3H)1/3 


(A.7) 


The  right-hand  side  (RHS)  of  the  vorticity  equation  consists  of 
stress-generated  expressions,  so  it  Bust  be  analyzed  more  carefully 
than  the  left,  with  particular  attention  paid  to  the  order  and  reference 
frame  in  which  time  Integrations  and  spatial  differentiations  are  per¬ 
formed.  This  is  crucial  for  proper  derivation  of  terms  of  the  form 


-r  x  -r  *  j(r»T)/P 
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-  I r  *  yr  •  |  )j(Htt* )dT* 

T 

*  j  VH  x  VH  •  [N(T-T,)Y(H,T,)]dT’ 
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where  j  and  y  are  stress  tensor  and  strain  rate  tensor,  respectively, 
and  the  subscripts  r  and  H  are  convenient  shorthands  to  specify  the 
particular  reference  frame  in  which  the  del-operator  (7)  is  defined. 
Begin  with  expressions  for  the  rate-of-strain  tensor  y  ^ (r,t) 
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despite  some  cumbersome  symbolism,  these  time-integrals  ere  evaluated 
for  constant  H,  while  with  in  the  stress  expressions  they  are  operated 


upon  In  the  (r,t)  eye tea-  Algebraic  manipulations  result  In  the  com¬ 
ponents  of  the  right-hand  side  of  the  vortlcity  equation 


7  l  (9  •  t)  (A.  12) 
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where  the  final  expression  comes  simply  from  direct  substitution  of  the 

definitions  of  f , ,  f„ ,  f_  and  £..  From  the  definitions  of  B  ,  X  and 
i  i  j  h  n  n 

Zr,  it  is  possible  to  simplify  these  expressions,  using 


“  nr*1”1  X  -  (n+l)r_(n+2)  Z  +  T 


(A.  13) 


32B  3T 

3^2-  -  n(n-l)rn~2  Xn  +  (n+1) (n+2)r"(n+3)Zn  +  (A.14) 


and  combining  like  terms  and  integrals,  the  stress-related  expression 


becomes  (A. 15) 
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The  differentiation  in  H  can  be  performed  inside  the  time  inte¬ 
grals  and  the  operations  applied  in  deriving  (A. 13)  and  (A. 14)  can  be 
re-applied  until  the  RHS  is  derived  in  the  form  which  was  employed  in 
the  final  T-equation  (A. 16) 
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Thus,  the  dimensional  vorticity  equation,  expressed  in  terms  of  the 
toroidal  field  function  T  (H,t)  can  be  written  as 


ST  T  2 

ar  -  2*2*  IT  -  vLl2  1  +  y  -  n(n+1>  pl 

+  r  {r  fW* 
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where  the  integrande  lj  were  given  previously  in  (111.41). 
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A.  2  Amplitude  Equations 


In  order  to  derive  the  boundary  conditions  at  the  bubble  surface* 
the  asymmetric  pressure  field  must  be  evaluated  ar  r  »  R,  (H-0)  as  must 
the  normal  and  tangential  extra  stresses. 

The  contribution  to  the  pressure  from  the  v-velodty  is  evaluated 
by  Integrating  the  equation  of  motion.  In  Lagrangian  coordinates*  for 
the  6-component 
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(A.  17) 


where  v„  *  (v  )_  .  Calculations  performed  in  the  development  of  the 

t7  V  0 

T-equation  give: 


(V  •  ;)Q  -  P  {  -p  H»(n+l)  +  1]  (A.  18) 
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The  angular  dependence  of  each  term  in  this  6-equatlon  except  the  pres¬ 
sure  term  has  been  determined  to  be  3T^/38.  Thus  the  pressure  term 
also  has  this  dependsncs  and  can  be  written  to  define  P_  as 


Pv(H,e,*,T)  -  Pv(H,t)  TfJ(e.d) 


Thus,  angular  dependencies  can  be  factored  from  equation  (A.  17)  and  an 
expression  for  Pv(H,t)  results.  Fortunately,  only  the  value  at  the 
bubble  surface,  H-0,  needs  to  be  explicitly  considered.  If,  In  the 
development  of  this  expression,  the  shorthand  notation 
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is  adapted,  then 
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(A.  21) 


The  stress  terms  are  also  considerably  simpler  When  evaluated  at  the 


surface  (A.  22) 
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Combining  these  results  (III. 44) 
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The  potential  flow  contributes  to  the  total  pressure  field  through  p^ 
which  was  derived  by  Plesset  to  be 


pp(r  -  R,6,4,t)  ■^(R»+3E»  +  2R  a)*JJ 


(A.  24) 


at  H  ■  0. 


Since  this  has  the  same  angular  dependence  as  pv>  an  analogue  to  Py  is 


found  to  be 
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The  unit  normal  to  the  surface.  g»  aleo  enters  into  the  0(e)  boundary 
conditions 


where  e  .  efl  and  e. 

*•  I  *"0  *  9 

tlons  are 
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are  the  unit  vectors.  The  dynamic  boundary  condi- 
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where  the  former  condition  is  for  zero  tangential  stress  and  the  latter 
expresses  the  balance  which  must  hold  between  any  discontinuity  in 
normal  stress  and  surface  tension  forces.  The  total  stress  tensor  o, 
includes  the  Isotropic  pressure  just  calculated  in  pp  and  pv  and  the 
extra  stress  t.  The  three  components  of  the  tangential  stress  condition 
are 

r  :  0  -  0 
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The  potential  flow  contributes  to  the  extra  stress  in  the  fluid  through 
the  rate-of-e train  tensor  for  that  flow  field 


(A.  28) 
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where  the  tine  dependent  function  Cn(t)  is  defined  ss 
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The  angular  tarns  in  curly  brackets  {  }  can  be  slnplifiad  using  the 
identity  for  spherical  harmonics 
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Thus,  for  the  $$-term 
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and,  for  the  6^-term 
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The  only  contributions  to  the  tangential  stress  condition  at  first 
order  in  e  will  cone  from  o^  or  c0r  at  0(ex)  and  orr,  o ^  or  o00  at 
0(e°) 
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These  combine  via  (A.26)  to  yield  the  tangential  stress  condition. 
Equation  (Ill . 46) . 
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The  normal  force  condition  at  first  order  in  c  reduces  to  only 
one  stress  component  balancing  the  0(e)  surface  tension  force,  but 
each  must  be  evaluated  at  the  deformed  surface 
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(A. 32) 
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The  stress  component  is  given  by 
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The  terms  needed  at  the  undeformed  surface  are  those  that  are  first 
order  in  c.  (A. 34) 
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The  first  term  can  be  replaced  by  an  expression  derived  from  the  tan¬ 
gential  condition  (III. 46). 

The  contribution  to  the  stress  resulting  from  the  surface  deforma¬ 
tion  is 
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at  H  •  0,  where  pQ  is  the  pressure  field  from  the  spherical  flow  given 
by 
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Contribution  (A.  35)  is  expressed  in  terns  of  a^,  R  and  derivatives  as 
an(T){12  pL2 

•CD 

-  pR(t)}  (e,4>)  <A*35’> 

Combining  (A.32),  (A. 34)  and  (A.35)  the  equation  of  notion  for  an(r) 
results  (111.45) 
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A. 3  Viscous  Limit 

A  necessary,  though  not  sufficient,  condition  for  the  validity 
of  these  results,  (111.41),  (111.45)  and  (111.46)  is  their  reduction  to 
the  expressions  for  a  purely  viscous  fluid  derived  by  Prosperetti 
(111.27),  (111.32)  and  (111.34).  This  Newtonian  limit  should  obtain 
for  the  special  choice  of  memory  function 

N(t)  ■  yfi(t)  Newtonian  (A. 36) 

where  6(t)  is  the  Dirac  delta  function.  This  relaxation  modulus  will 
serve  to  reduce  the  time  integrals  to  the  value  of  the  integrand  at  the 
time  when  the  argument  of  the  delta  function  is  zero. 

[  u6(T-T,)f(T,)dT’  »  yf(x)  (A. 37) 


Thus,  the  distinction  between  Eulerian  and  Lagranglan  time  Integrals 
vanishes  since  the  variation  of  geometric  parameters  such  as  R(t)  over 
past  times  is  no  longer  relevant. 

The  simplifications  that  the  Newtonian  stress  relation  allow  in 
the  model  system  are  substantial.  For  the  T-equatlon  in  the  form  of 
(111.41),  the  elastic  integrals,  those  multlpled  by  Go,  are  identically 
zero.  The  cancellation  of  terms  in  the  Integrands  Ip  Ij,  and  which 
occurs  to  leave  only  the  term  multiplied  by  v  in  (III. 41)  is  best  seen 
in  expression  (A.16).  The  distinction  between  terms  like 


|TN(T-T,)Xn(H,T,)L,<n“2)dT' 


and 


L  J  N(T-T,)Xn(H,T,)L,(n“5)dT' 


i 

i 
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j 

■3 
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1 
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disappears  for  the  purely  viscous  case  since  both  are  now  equal  to 
vXn(H,T)L(n"4) 

Terms  containing  and  ZR  cancel  completely,  and  only  the  terms  in¬ 
volving  Tn  remain.  It  is  then  straightforward  to  transform  the  Lagrangian 
purely  viscous  T-equation  to  the  Eulerlan  result  (III. 27). 


f  +  I? 


W0  -  »<"«> 


The  tangential  stress  condition  (111.46)  is  even  more  easily  re¬ 
duced  to  the  desired  viscous  limit.  The  distinction  between  terms  such 
as 


j  f  "•  jr 


and 


r  .*»• 

N'  dr’ 


vanishes  in  this  case  and  (III. 32)  obtains  directly 


The  purely  viscous  no real  stress  balance  is  more  difficult  to 


produce  from  the  non-Newtonian  expression.  First*  Fy  (111. 44)  Bust  be 
obtained  for  the  Newtonian  Modulus  in  Eulerian  coordinates*  then  the 
result  can  be  employed  in  the  amplitude  equation  (III. 45).  In  equation 
(III. 44)  most  of  the  tens  in  the  two  tlae  Integrals  will  cancel  in  the 
viscous  Halt.  When  there  is  no  externally  istposed  flow*  the  result 
is 

viscous 

Pv(H-0,t)/p  -  Pv(v-R,t)/p  (A. 38) 


.  _v  II  +  Pi**!?  R“ 

V  3r  n+1  R 
00 

+  R  Tn(R,t)  -  Rn  j  S 
R 


<R)or 


~n  r3T  a4U2  3TtHc 
ta?  +  R(s)  ?slds 


Using  the  viscous  T-equatlon  (III. 27) 

If  "  '  fe  ^<f>2T]  +  -»(■*>  f*> 

and  spatially  Integrating  by  parts 
Py(r-R,t)/p 

•e 

■  n(f)  I  I<f>! '  »)<§>“  «•.*)« 

which  is  equivalent  to  (III. 35).  With  this  result  the  viscous  equation 
(111.34)  emerges  from  (III. 46). 


APPENDIX  B  :  0(e2)  Terms 
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B.l  Introduction 

The  validity  of  the  present  analysis,  which  is  restricted  to 
"small  e"  obviously  depends  upon  that  small  parameter.  One  way  to 
evaluate  the  radius  of  convergence  of  the  enalysis  is  through  experi¬ 
mental  trials  and  these  seem  to  indicate  good  agreement  for  deforma¬ 
tions  which  are  substantial  fractions  of  the  overall  bubble  radius 
(see  Figures  24-27).  Previously  published  comparisons  of  the  non¬ 
linear  analysis  of  Chapman  and  Plesset  (1971)  with  linearised  results, 
these  for  invisdd  fluids,  are  encouraging  since  they  show  good  agree¬ 
ment  until  the  final  stages  of  collapse.  Analysis  can  also  give  some 
estimate  of  the  range  of  validity  of  a  linearisation  without  the  need 
for  solution  of  the  full,  nonlinear,  problem  through  the  generation  of 
the  neglected  terms  and  determining  the  conditions  under  which  they  are 
small  compared  to  the  terms  which  are  retained. 

The  nonlinear  terms  which  are  dropped  in  the  derivation  of  the 
T-equatlon  (III. 41)  from  the  exact  vortldty  equation  (22)  are  terms 
of  second  order  in  e.  Such  second  order  terms  are  also  quadratic  in 
the  spherical  harmonics,  which  removes  the  separability  of  the  equa¬ 
tion  into  an  (r,t)-dependent  function  and  an  angularly-dependent  one. 
This,  in  turn,  prevents  the  manipulation  which  totally  removed  the 
angular  dependence  from  the  final  expression  (see  equations  A. 5  &  A. 13) 
and  the  resulting  independence  of  this  equation  from  those  for  any 
other  value  of  the  index  n.  Thus,  the  explicit  evaluation  of  the 
next  term  in  the  expansions  for  velocity  or  pressure  (equstlons  III. 16, 
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17)  Is  considerably  aore  complicated  than  that  for  the  0(ex)  term. 

Since  both  the  full  nonlinear  problem  and  the  0(e2)  problem  are  so  dif¬ 
ficult  ,  the  error  introduced  by  the  nonlinear  terms  is  estimated  by 
evaluating  them  for  the  values  obtained  from  the  linear  analysis. 


B.2  Convection  and  Inertia 

The  nonlinear  terms  in  the  LHS  of  the  vortidty  equation  which 
were  neglected  are 


<YV  +  Yp>  <yv  +  v  >  (b.i) 

and  represent  convected,  inertial  phenomena.  They  also  bring  the  poten¬ 
tial  velocity  field  Vp  into  the  analysis  for  vortidty.  The  expressions 
for  each  component  of  the  velocity  and  vortidty  are  given  in  Section 
III  and  again  in  Appendix  A  (III. 19,  III. 28,  A.l,  III. 21,  A.4).  A 
double  summation  over  the  lower  index  on  the  spherical  harmodcs  is 
now  necessary.  A  new  symbolism  is  adopted  to  insure  the  explicit  reten¬ 
tion  of  both  values.  Let  "k"  be  one  subscript  and  "1"  the  second  and 
let  v  be  the  total  0(c)  velocity. 

Y  -  Yv  +  Yp  (B.2) 

Then,  this  means 

<v.  •  «t>k+2  <\ + v  !>  ♦ 


(B.3) 


For  convenience »  the  potential  contribution  is  represented  through 


(^(r.t)  -  (*)k+2  (ijj  +  ^  j[>  <B**> 

This  symbolism  can  be  used  to  express  the  0(e2)  terms  in  the  Eulerian 
frame  as 


»l>r 


(B.5) 
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(B.6) 


lYk  *  v“t3e 
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(B.9) 


8Bk  8  Ti  8Yi 
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,»T*»V  i  3\  *\ 
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(B.10) 
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■i 


T.  -  8Y.  B(  C.  .  .  8Y„ 

r_£]  {-  — —  — — r  (_  4.  _JL)  1 - L_  ~i 

lr*Jl  sin8  8*  lvr  l+l 1  86  sin6  8*  J 


8Y.  Bf  Cf  8*Yf  8Y 

+  ’55"  1  +  t*T)  ’■W  *  “te  W* 


-  «.--5T+ W» 


A  tractable  expression,  and  one  which  ia  etlll  Meaningful .  re¬ 
sults  If  axlsywetry  is  assumed  and  only  one  node  ia  present*  a.g. 
k  ■  1*2.  The  r-  and  8-cooponents  of  the  cross  terms  are  identically 
aero  and  the  *-tera  is  (B.9»  10) 


3BV  1  3TU  Tu  SYU 

“  <*k-  3?  +  W-7aT+^  ar 

B.  C.  X.  3Y.  32Y. 

+  j.  _£  \  (  *\  _ k _ JL 

'r  k+l' V"  36  36s 
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3Y, 


A  ,  Sc .  ro.  ^k 

(r  ){r  +  k+15  36  cote  36 


SB,  T.  3Y. 
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3B.  2T.  ,  3T.  3Y. 

(Tk  "  3T  +  Sc)(“^  "  r  WL)Yk  3T 


B,  CL  T.  3ZY  3Y 

+  (T  +  th0  <7?)  (30^  -  “te  -jg-y 


»V?\ 

36 


From  the  sample  calculation  of  Appendix  C.2  values  for  each  of 
these  functions  can  be  assigned  and  the  relative  magnitude  of  the  convec¬ 
tion  terms  evaluated,  e.g.  at  r  ■  K  (B.ll) 

Step  #10:  time  ■  2.50  x  10”s  sec 

R  -  2.42  x  10"1  cm 
R  ■  1.47  x  10*  cm/sac 
a  ■  9.5  x  10-2  cm 


a  -  -1.863  x  10*  cm/sec 
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a2  ■  -8.70  x  10* 

|^(R)  -  -4  x  107  sec”1 
3r 

T(R)  ■  -1.14  x  10*  cm/sec 


Step  #11:  time  ■  2.75  x  10“s  sec 

R  -  2.46  x  10"1  cm 

R  ■  1.41  x  ios  cm/sec 

a  ■  9.1  x  10“2  cm 

a  -  -1.943  x  102  cm/sec 
a2  -  -8.805  x  10* 

|^(R)  -  -1  x  10*  sec"1 
T(R)  -  -1.20  x  10*  cm/sec 


Also  at  r  -  R 


3B2 

T - -  a  0 

l2  3r 


C2  ■  a  +  2a  *v  10*  cm/sec 


r  •  r Rn_1  °2 ' 4  x  10> 


Thus  the  factors  in  the  second  order  expression  ere 


32 

(T2  -  -jj—  +  C2)  j  ^  10s  cm/sec 


r-R 


T  3x 

<2  ~  -  |  ^  2  x  10*  cm_1sec-1 
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B  C 

(~-  +  ^~)  |  a.  4  x  10s  cm/sec 
r-R 


(B.12) 


T2 

(jt) 


2  x  io%  sec“Icn“1 
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3Y, 
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(B.13) 


3%  8Y  3Y« 

+  8  x  107(-^  -  cote  jf)  jf  )sec“2 


These  second  tens  ere  to  be  coapered  with  the  inertiel  tens  et 


first  order  in  t  Which  ere 


1 


.  [i  U  (10  +  I  (5  SI  .  2E  I>1  ^ 

lR  3t  w  SlKJr  R;  J  36 


(B.14) 


Numerical  data  gives  this  expression  (B.14)  a  value  of  about  |3  *  10 1 1 | 
3Y2/30  at  r  ■  R.  For  these  terms ,  at  first  order  in  e,  to  dominate  the 
nonlinear,  second  order  terms  of  (B.13) 


e("B.14")  >  e2  (HB.13'') 


(B.15) 


or  numerically 


1 3  *  1012l 

"  2  x  10ll| 


(B.16) 


To  arrive  at  this  result,  the  common  factor  dY^/ae  has  been  cancelled 
and  the  other  angular  terms  in  (B.13)  neglected  since  they  are  0(10°) 
or  less;  they  consist  of  combinations  of  sine  and  cosine  functions. 

For  large  R  values,  as  in  this  case,  these  inertial  terms  are  dominated 
by  the  contributions  of  and  3T/3r .  When  this  occurs  the  requirement 
on  c  becomes 


« (rr  +  ~>  *  1 

2a  R 


(B.17) 
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B.3  Mon-Newtonian  Nonlinaarltiea 

A  second  instance  for  which  the  importance  of  second  order  terns 
must  be  evaluated  does  not  arise  from  a  straightforward  expansion  of 
existing  tens  in  the  vorticity  equation,  but  from  a  particular  choice 
of  constitutive  model  for  the  stress.  Cross  tens  can  arise  from  the 
nonlinearities  of  convection,  co-defomatlonal  or  co-rotatlonal  tens. 

One  method  by  which  to  approximate  these  contributions  while 
maintaining  the  fon  of  the  previous  results  is  to  employ  quasi- linear 
expressions,  in  the  sense  of  Bird,  Armstrong  and  Hassager  (Chapter  7). 

A  Taylor  series  expansion  in  convected,  co-rotational  or  co-deformational 
coordinates  can  be  employed  to  approximate  the  difference  in  stress 
each  formalism  would  predict.  In  the  integral  constitutive  model  the 
strain  rate  at  time  x',  j(x)»  or  its  equivalent  can  be  replaced  by 


I 


co.vect.d(T,)  '  *(t>  +  jtt  <t> 

-  Y<x)  +  (t-t*)  (|^  +  ey  *  VY)| 


T«X 


(B.18) 


in  the  convected  case.  The  present  0(cl)  model  already  incorporates 
an  evaluation  scheme  which  is  better  than  the  first  two  terms  in  the 
second  expression  of  (B.18)  since  the  actual  value  of  strain  rate  at 
past  times  is  evaluated  and  Included  in  the  analysis.  If  this  evalua¬ 
tion  is  termed  v _ »(?')  then  another  approximation  to  the  con- 

vected  strain  rate  at  past  times  is 


l-ati. 


H-ll 


^convected^1  ^  "  ^present  ^ 


+  e  (t-t ’ )V  •  Vy 


(B.19) 


Similar  expressions  arise  when  co-rotational  contributions  are  included 


(B.20) 


ico-rot(T,)  "  ?(t)  +  e(T“T’>fei(T) 


?convected^T  ^ 


+  c(t-t')  J  K«  *  j)  ”  (y  *  “) 


for  the  vorticity  tensor  u.  And  also  when  co-deformational  contribu¬ 
tions  are  considered 


?co-def(T,) 


(B.21) 


Y  (T«\ 

•convected' 


+  e(t-T')I(Vv)T  •  y  +  j  •  (Vy) ] 


These  expensions  are  linear  in  time  and  can  only  be  accurate  as  long 
as  the  time  derivatives  do  not  change  sign. 

All  of  these  instantaneous  strain  rates  will  be  Integrated  over 
time  in  the  same  stress  expression.  Thus,  if  the  instantaneous  dif¬ 
ference  between  the  present  strain  expression  and  those  including 


■  v-T  ■:  . 


convected,  co-rotational  and  co-de format ional  contributions  is  small, 
the  resulting  stress  differences  will  also  be  small.  This  contrasts 
sharply  with  any  comparison  between  purely  viscous  and  viscoelastic 
stresses  because  they  cannot  be  validly  approximated  by  taking  only 
taneous  strain  rates  into  account. 

By  reverting  to  the  Taylor  series  expansion  for  the  strain  rate 
as  presently  employed,  three  terms  emerge  which  must  be  less  than  the 
present  strain  rate 


iconvected*  ^co-rotational* 


and  Ay 

ieo-i 


* co-deformat ional 


(B.22) 


where 


<  ^present  "  ?<’>  +  ^ 


4Xconvected  *  *<"'>?  '  'X 


Aico-rot  “  I 


Aico-def  "  e<T"T'>t(?v)T  *  y  +  J  •  <?y>J 


The  convected  difference  terms  can  be  represented  by  the  (rr) 


component 
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<y  •  ?j>rr  -  yt  a?  i 
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v 
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+  r”  36  Irr  +  rsin6  ?rr 
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(B.23) 


Combining  this  expression  with  those  for  the  full  velocity  vector 
and  rate  of  strain  tensor  of  Appendix  A,  the  derivatives  are  all  on  the 
order  of  the  particular  strain  component  divided  by  the  local  radius , 


e-g* 


—  (y  )  -  2F—  (T  -  — ) 
3r  «rr'  ^ar*'1  3r' 


<*>k+2  V  \ 


0.24) 


emerges  from  (A. 9)  and  (A. 28).  It  can  also  be  shown  from  the  definition 

of  B  (A. 2)  and  (A. 13)  and  (A.14)  that 
n 

|f2(T  -  |f)  -  -«k-l)Ck-2)Vk‘3 

+  (k+l)(k+2)0t+3)Zkr'<kM)  0.25) 

♦  4  isth*tn 

**  2k+l 


which  means  that 


(B.26) 


is 


In  the  evaluation  of  the  present  strain  rate  the  (rr)  component 


,3T  92Bv  _m 
'3r  "  3r*;  n 


(A.  9) 


-  -2[n(n-l)rn”2X  +  (n+1)  (trt-2)r“(n+2)Zl 
n  n 


from  (A.  14).  The  sample  calculation  assigns  numerical  values  to  the 
terms  in  (B.22),  (B.23),  (B.26)  and  (A.9),  at  r  -  R  (see  B.ll) 


|j|  i  1.75  x  10s  sec**1 

3-r 

*  8.4  x  io#  sec"2 


R  *  2.4  X  10"1  cm 


(B.27) 


1  v J  *  10s  sec/sec 


For  this  case,  the  convected  contribution  to  the  strain  rata  is  smsll 


(B.28) 


C<T— T'KIO’)  jilVlff! 

1.75  «  10 5  +  «  10*  )  *  1 

which  is  true  for  all  t*  <  t  whan 

e  <  1.15  (B.281) 
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A  similar  analysis  can  be  carried  out  for  the  co-rotational  terms. 
For  the  axlsymmetric  case  the  necessary  calculations  show  that 


(B.29) 
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has  only  four  non-sero  te 


,  which  are  given  by 
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The  expressions  for  the  vorticity  tensor  (A.4)  show  thet 


to 
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(B.28) 


And  thus  the  e  condition  for  co-rotetion  becomes 

e (t-t  1 ) J~|  |i| 

- r-  <  1  (B.29) 

i  +  (T-f)  ^ 

3t 

which  numerically  Implies,  for  this  case 

t  <  0.96  (B.291) 

Finally,  for  the  co-deformatlonal  formalism  the  only  complicstion 
is  the  relevance  of  both  the  potential  and  v-velodty  fields.  The 
velocity  gradient  consists  of  components  of  the  form 

W  *  0  (|2  +  h  (B.30) 


and  the  e -condition  emerges  as 
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e  <  0.8  (8.311) 


B. A  Amplitude  Equations 

The  Inaccuracy  of  the  linear  approximation  applies  to  the  non- 
Newtonian  expressions  in  the  boundary  conditions  just  as  it  did  for 
the  vortlcity  equation.  The  results  are  the  same  as  for  T,  since  the 
same  stress  expressions  are  used  to  evaluate  these  quantities  at  the 
surface.  However  there  Is  one  place  where  these  0(e2)  terms  *  ould  become 
Important  for  the  boundary  conditions  in  addition  to  those  for  the 
T-equation.  The  unit  normal  to  the  surface  used  to  formulate  these 
conditions  Is 

_  r  .  f  jW.  ,  *n  3  . 

5  R  sr  e  r  38  *e  e  sine  3+ 


The  correction  to  this  expression  at  second  order  In  c  Is 


-f  +‘ 


1  3Yv2 

sine  3*' 


(B.32) 


Thu*  the  correction  will  be  eaall  if 


e 


< 


2 


R 


or 


e  <  5 


for  these  numerical  values. 


APPENDIX  C 


C.l  Method  of  Solution 


The  system  of  equations  governing  the  toroidal  field  To(H,t) 
and  the  nonapherical  amplitude  a^  can  be  solved  through  a  series  of 
relatively  straightforward  finite  difference  calculations  combined 
with  linear  algebraic  techniques.  The  basic  procedure,  which  ellows 
simultaneous  solution  without  iteration  begins  with  the  calculation  of 
«  "particular''  and  a  "homogeneous''  solution  to  the  vortlcity  equation 
for  T  .  The  correct  linear  combination  of  these  solutions  is  deter- 

Q 

mined  using  both  the  amplitude  equations  and  finite  difference  approx¬ 
imations  to  the  nonspherlcal  velocity  and  acceleration.  A  linear  sys¬ 
tem  of  three  equations  is  solved  which  determine  amplitude,  velocity 
and  the  proper  linear  combination  simultaneously.  Internal  consistency, 
l.e.  convergence,  of  the  solution  was  checked  by  systematic  variation 
of  numerical  parameters  and  algorithms  within  each  program  segment. 

The  symbolism  employed  reflects  the  restrictions  of  computer 
output  and  also  the  discretisation  of  continuous  functions  performed 
in  applying  finite  difference  approximations.  Spatial  variation  is 
signified  by  the  index  "1"  and  temporal  dependence  by  a  second  index 
"J".  The  spatial  grid  was  generated  by  specifying  the  Lagranglan 
coordinate  H(I)  as  a  geometric  series.  An  initial  H  value  was  speci¬ 
fied  to  be  HI  and  a  ratio  ph  was  also  chosen  s.t. 


•  P 


H 


(C.l) 


vhere 


H(l)  -  HI 


•nd  H  -  0 

o 

The  ratio  was  leas  than  two  (2.0)>  since  this  was  found  to  be  neces¬ 
sary  for  good  approximations  of  spatial  derivatives,  and  greater  than 
unity.  This  lower  bound  served  to  concentrate  the  spatial  points  in 
the  vicinity  of  the  bubble  surface  while  still  allowing  relatively  few 
steps  to  span  a  large  spatial  range.  This  large  range  is  desirable  for 
the  approximation  of  the  integral  term  a&  (see  Equation  III. 30) 


aQ  -  |  S_n  T(S,t)dS 
R 

r 

max 

i  |  S‘n  T(S,t)dS 
R 


(C.2) 


Trial  and  error  calculations  demonstrated  that,  while  the  value  of 
r  necessary  for  convergence  is  not  large,  it  would  require  hundreds 
of  equally  spaced  steps  of  HI  to  reach  that  value  when  Hi  is  small 
enough  to  insure  good  approximation  of  derivatives  near  H-0.  Thus, 
this  geometric  series  was  used.  Convergence  could  usually  be  achieved 
for  values  such  as 


R  -  10-*  cm 
HI  -  10-11  cm* 

(C.3) 

PH  -  1-5 
%  -  25 

where  Hg  is  the  total  number  of  spatial  steps.  For  these  parameter 
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values  r _ *  10” 2  ca. 

The  sero-order  problem  of  determining  the  radius  profile  R(t) 
was  solved  In  a  straightforward  manner,  employing  standard  predictor- 
corrector  methods.  The  evaluetion  of  the  elastic  contribution  to 
stress  was  simplified  by  the  recursive  property  of  an  exponential 
memory  function  explained  below.  A  Simpson's  rule  approximation  of 
the  stress  integral  results  in  (C.4) 

|  exp  ~  £  -  AU.)  4rr(t‘)dt' 

—00 

.  ..  *rr(t  +  4t>  +  exl>('  r>*rr(t) 

*  “ - 2 - 

t 

+  exp(-  |  ^(f)  exp(^)dt' 

— «0 

e  •• 

The  R  profiles  so-generated  were  then  available,  along  with  R  and  R 
values ,  as  inputs  to  the  0(c)  solution.  Alternatively,  the  profile 
generated  by  the  data  fit  explained  in  Appendix  D  was  also  available. 

As  in  the  numerical  solution  of  many  P.D.E.  when  time  is  among 
the  independent  variables,  the  procedure  was  to  evaluate  or  specify 
all  quantltltes,  e.g.  XTn(X,J)  A  To(H,t) ,  for  I  »  1,  2,  ...,  Kg,  at 
some  time  step  "  J" •  then  use  those  values  to  generate  the  corresponding 
numbers  for  the  next  time,  "J  +  l".  The  finite  difference  approxima¬ 
tion  to  the  T-equation  employs  a  central  difference  in  the  time  deriva¬ 
tive  and  a  mixed  implicit-explicit  weighting  of  the  remainder  of  the 
equation,  f^,  at  J  and  J  +  1.  The  function  f^  results  when  the 


vorticity  aquation  la  expressed  In  the  for*  (C.5) 


(C.5) 

i  £t(H;  R,R,T,  II. 12 .13) 

The  weighting  was  characterized  by  a  relaxation  parameter 
0  <_  1,  where  $T  ■  1  corresponds  to  a  completely  explicit  algorithm. 

This  relaxation  parameter  was  among  the  numerical  parameters  varied  in 
early  calculations  to  insure,  and  then  optimize,  convergence.  Spatial 
derivatives  were  approximated  by  a  two-point  central  first  derivative 
and  three-point  second  derivative.  These  were  chosen,  despite  the 
resulting  requirement  of  small  spatial  steps,  because  the  resulting 
matrix  form  of  the  XT-equatlon  allows  a  non-iterative  solution  by 
Thomas'  method. 

The  algorithm  generated  from  these  choices  was  applied  twice. 

A  "particular"  solution,  TZ(I)  at  (J  +  1),  results  when  XT(I,J)  is  re¬ 
tained.  A  second  "homogeneous"  solution,  DTF(I)  at  (J  +  1),  is  pro¬ 
duced  when  the  previous  XT  values  are  neglected.  The  linearity  of 
the  equation,  allows  the  general  solution  st  (J  -f  1)  to  be  expressed 
es 

TT(I,  XRAT)  -  TZ(I)  +  XRAI  *  DTP(I)  (C.6) 

where  the  fsctor  XRAT  must  be  determined  by  boundary  conditions.  The 
"homogeneous"  and  "particular"  functions  ere  also  sufficient  to 
determine  the  integral  on  and  pressure  term  Py  within  this  same  unknown 


factor  XRAT.  The  no— nclatura  employed  la  defined  by 


<-  an  4  UI1  -  UI1Z  +  XRAT  *  Ml 

and  p 

Pv(H-0,t)  •  PRZ  +  XRAT  *  DPR 

where  U11Z  and  PRZ  obtain  for  the  particular  aolution  TZ  and 

Dll  and  DPR  from  DTF. 

The  tangential  atress  equation  (III. 46)  can  be  recaat  as 

y  (J  +  1)  -  f .(a  .  TZ,  DTF,  D11Z,  Dll;  XRAT;  R,R; 

XI  A  XI 

(C.7) 

past  values,  physical  parameters) 

where  all  the  independent  variables  preceding  the  first  semi-colon  are 
for  time  (J  +  1)  and 

yn(J  +  1)  -  an(J  +  1)  (C.8) 

The  normal  force  condition  is  somewhat  more  complicated,  yielding 

yn(J  ♦  1)  -  in(J  +  1)  (C.9) 

•  fy(fn»  *n»  12  *  DTF, 

PRZ,  DPR;  XRAT;  R,  R,  R; 

past  values,  physical  parameters) 
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Various  dosed  integration  algorithms  are  available  to  give  finite 
difference  approximations  to  these  equations.  Here*  four  different 
algorithms  were  employed  —  explicit;  simple  two-point  mixed  implicit 
explicit  (eA*By);  mixed  two  and  three-point  (B^,  By»  6^,  B^)  and  a 
modification  of  a  Hamming  predictor-corrector  employing  four  points. 
These  were  all  tried*  and  choices  made*  once  again  to  insure  and 
optimize  convergence.  All  algorithms  lead  to  linear  equations  in 

y  (J  +  1),  a  (J  +  1)  and  XRAT 
n  n 


yn  *  WA(L)  +  an  *  WB(L)  +  XRAT  *  WC(L)  -  WD(L)  (C.10) 

j 

i 

A 

for  L  -  1,  2  and  3,  and  where  the  coefficients,  HA,  WB.  HC  and  HD*  de¬ 
pend  only  on  past  values*  physical  parameters  and  R*  R,  ft,  TZ,  DTF, 

UI1Z,  Dll,  PRZ  and  DPR*  l.e.  they  are  effectively  constants  with  respect 
to  Y^,  aQ  and  XRAT.  This  is  a  well-posed,  determinate  system  which  is 

then  solved  for  a  (J  +  1),  Y  (J  +  1}  and  XRAT* 
n  n 

With  these  values,  all  quantities  can  be  updated  and  the  process 
repeated  for  succeeding  times. 

C.2  Sample  Calculation 

This  method  of  solution  is  bast  Illustrated  by  a  sample  calcula¬ 
tion.  For  a  purely  Newtonian  fluid  of  viscosity  0.10  cp,  and  the  values: 


time  step  -  2.50  x  10  sec 


HI  -  10"*  cm* 
pH  -  1.50 
%  “  30 

initial  radius ,  R  -  1.98  *  10" 1  cm 
o 

initial  amplitude!  aQ  *  1.27  x  10_1  cm 
initial  nonspherical  velocity,  aQ  ■  -4.0  x  102  cm/sec 
The  values  at  time  step  #10  were  found  to  be 
time  -  2.50  x  10~5  sec 
R  «  2.42  x  10“l  cm 
R  -1.47  x  10*  cm/sec 
R  -  -2.21  x  io7  cm/sec2 
XT(R-O)  -  -1.140  x  10s 


I  - 

1 

2 

3 

4 

-1.140 

-1.141 

-1.142 

-1.144 

w 

10 

15 

20 

25 

H 

o 

N 

. 

1 

-1.712 

-2.426 

-2.105 

a  «  -8.702  x  10s 
a  -  9.553  x  10-2 
a  -  -1.863  x  10s 

To  calculate  these  values  at  time  steo  #11,  first  these  values  are 
generated. 

t  *  2.75  x  10“5  sec 
R  -  2.46  x  10-1  cm 
R  ■  1.41  x  io3  cm/sec 
R  •  -2.19  x  107  cm/sec2 


and 

TRZ  -  -5.315 

DTRF  •  0.996 

U11Z  - 

-8.804  x  10s 

Dll  -  6.422  x 

* 

1 

o 

rt 

0I2Z  - 

-2.933  x  10® 

DI2  »  2.804  x 

102 

so  that 


WA  WB  WC  WD 

L  -  1  3.39  x  10'2  -6.42  x  io1  -1.27  x  10“2  -5.71  x  io1 

2  1.016  x  io6  8.35  x  101  1.62  x  io2  - 

3  9.37  x  IQ"7  -1.00  x  10°  -  -9.26  x  10~2 


which  can  be  used  to  determine 

a  -  -1.948  x  10s 
a  -  9.073  x  lo-2 
XRAT  -  -1.17  x  10s 
and  XT  (H-0)  -  -1.170  x  10s 


I  - 

1 

2 

3 

4 

5 

-1.171 

-1.172 

-1.173 

-1.175 

-1.177 

XT(I) 

10  3 

10 

15 

20 

25 

-1.234 

-1.732 

-2.499 

-2.162 

Variation  of  time  step,  spatial  step  and  number  of  spatial  steps 


was  performed  with  convergence  to  these  amplitude  values 
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APPENDIX  D 

Data  Analysis 

D.l  Shape  and  Site  Fit 

A a  stated  in  Section  V,  from  each  photograph  the  equivalent 
radius  R,  and  some  measure  of  the  nonsphericity  of  the  bubble  is  de¬ 
sired.  The  shape  of  these  axisyometric  bubbles  is  assumed  to  be  well 
fit  by  a  series  containing  up  to  three  nonspherical  modes 


r(0,d,t)  -  R(t)  +  a2(t)  y£(0,*) 

+  a3(t)  Y°(0,*)  +  a4(t)  yJ(0,*) 


(D.l) 


where 


Y°(0,*)  -  P2(cos0)  -  7  (3  cos2  0-1) 

Y°(0,*)  -  ?3(cos0) 

-  j  (5  cos*©  -  30cos  ) 

Y®(0,4)  -  P4(co.0) 

-  i  (35  cos*6  -  30  cos20  +  3) 

In  order  to  fit  this  function  to  the  actual  bubble  image  it  is  also 
necessary  to  determine  the  best  location  for  the  origin,  and  center- 


J 

i  J 
d  ' 


i 

.1' 


D-2 

line  axis  from  which  r  and  8  can  be  measured. 

The  basic  strategy  by  which  a  fit  can  be  made  does  not  vary  with 
a  particular  procedure.  This  strategy  Involves  digitally  encoding  the 
location  of  points  on  the  surface  of  the  bubble  in  the  chosen  coordinate 
system,  then  performing  a  least  square  analysis  to  optimise  the  param¬ 
eter  values. 

1)  The  simplest  procedure  employs  just  two  measurements  of  bub¬ 
ble  image,  the  major  axis,  which  is  assumed  to  be  the  largest  horisontal 
dimension  and  the  minor  axis,  the  vertical  dimension.  These  occur  at 
6*0  and  ir  and  0  *  and  ,  respectively  when  the  contribution  of  the 
third  and  fourth  modes  is  assumed  to  be  neglibible.  If  the  lengths  of 
the  two  axes  are  labelled  £q  and  lv/2,  then  values  can  be  determined  for 
R(t)  and  a2(t)  through 


£o  -  2R  +  2a2  P2(l) 

«  2R  +  a2 

£ir/2  «  2R  +  2a2  P2<0) 
-  R  -  a2/2 


which  means  that 

i  +  2  tv/2 

n _ O _ 


‘o  . 

•2  •  r  *  * 

where  a^  *a^  -  0  has  been  assumed. 
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This  first  procedure  cen  be  accomplished  directly  free  the  photo¬ 
graph  or  eore  accurately  hy  magnifying  the  image.  Magnification  mas 
accomplished  by  opaque  projection.  The  aucceeding  three  procedures 
require  the  additional  accuracy  afforded  by  the  projection  technique. 

2)  By  guessing  the  centerpolnt  of  the  image  and  the  axis  of 
symmetry  a  polar  coordinate  digitisation  becomes  possible.  The  (r,6) 
coordinate  pairs  obtained  can  then  be  fit  by  a  linear  least  squares  fit 
to  the  function  (D.l). 

Unfortunately,  the  lack  of  fore-aft  symmetry  (right-left)  in 
photos  such  as  Figure  28(b)  and  28(f)  makes  the  a  priori  estimation  of 
the  origin  difficult.  The  location  of  that  point  also  needs  to  be  fit. 

Figure  28  also  contains  frames  which  display  another  feature  of 
the  Images  which  makes  analysis  difficult.  The  bubble  appears  to  have 
flattened  ends  in  frames  (e),  (f),  (g)  and  (h).  The  ends  are  not 
really  flat,  but  are  the  reault  of  the  indentations  Illustrated  in 
Figure  9(d)  and  9(e)  which  do  not  appear  in  the  photos  since  the  interior 
of  the  bubbles  is  not  visible.  Thus,  editing  of  the  digitised  data  may 
be  neceasary  for  aome  frames.  Unlike  most  fitting  procedures  wherein 
greater  confidence  can  be  placed  in  results  for  larger  numbers  of  points, 
here  judgement  of  "flat  spots"  is  required. 

The  two  methods  expected  to  be  moat  accurate  include  a  fit  to  the 
origin.  They  are: 

3)  The  centerline  corresponding  to  the  axis  of  symmetry  la  fixed 
a  priori ,  then  digitised  data  is  fit  to  generate  the  best  value  of  ft, 

*2  and  the  centerpolnt  on  the  axis. 

4)  Same  as  method  3),  with  the  addition  of  s  fit  to  a^  and  a^. 
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Both  procedures.  3)  and  4).  employed  the  sane  data.  Careful 
estimation  of  the  centerline  vaa  made  on  enlarged  Images,  which  was 
aade  possible  by  the  excellent  angular  (up-down)  syametry  of  the  out¬ 
lines.  A  first  guess  of  the  centerpoint  was  aleo  aade.  Froa  this 
center  lines  were  drawn  radially  at  15*  to  20°  Intervals.  At  the 
intersection  of  these  line  segments  with  the  surface  image,  rectangu¬ 
lar  coordinates  were  aeasured.  Rectangular  coordinates  were  eaployed 
because  this  geometry  allows  simple  transformations  when  the  center- 
point  is  translated  along  the  axis  to  Improve  the  data  fit.  The  best 
values  of  the  desired  amplitudes  were  then  determined  by  an  iterative 
nonlinear  least  squares  algorithm  which  minimised  the  sum  of  the  square 
residues.  Trials  were  run  with  all  the  data,  and  also  with  editing 
for  "flat  spots". 

Results  for  methods  1),  3)  and  4)  are  shown  in  Table  Dl.  Here 
results  are  shown,  not  only  for  different  fit  procedures,  but  also  for 
different  times,  different  photos  at  the  same  time,  and  for  data  with 
and  without  editing.  Looking  first  at  R  values,  the  variation  which 
occurs  between  methods  3)  and  4)  is  seen  to  be  minimal  for  different 
analyses  of  the  same  photo.  The  reproducibility  between  different 
photos  is  not  as  good,  with  a  total  variation  of  about  10%  for  the 
five  photographs  analysed  at  750  ysec  for  methods  3)  and  4).  Method  1) 
doea  not  agree  quite  as  well  with  the  other  two  methods  eaployed;  the 
R  values  from  method  1)  are  consistently  5-10%  smaller  than  the  others, 
but  do  show  proper  trends. 

The  *2  values  show  mors  variation  then  the  t  values  due  to  fitting 
procedure  and  to  photograph  to  photograph  reproducibility.  There  is  no 
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*2 
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% 
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OO-1!*) 
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2 

3.00 

0.67 

3 

U 

3.1* 

1.06 

4 

3.17 

0.96 

0.63 

-0.40 

3 

1* 

3.17 

1.24 

4 

3.1* 

0.66 

0.63 

•0.2* 

700 

a  1 

2 

3.25 

1.10 

3 

U 

3.3* 

1.21 

4 

3.3* 

0.95 

-0.79 

0.30 

b  1 

2 

3.20 

1.20 

3 

10 

3.33 

1.19 

4 

3.32 

1.07 

0.04 

0.29 

c  1 

2 

3.30 

0.76 

3 

IS 

3.51 

0.74 

4 

3.52 

0.65 

0.17 

0.17 
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a  1 

2 

3.3* 

0.33 

3 

IS 

3. SB 

0.36 

4 

3.58 

0.36 

-0.19 

0.01 

3 

17 

3. 56 

0.38 

4 

3.56 

0.37 

-0.22 

0.04 

3 

IS' 

3.56 

0.42 

4 

3.56 

0.38 

-0.22 

0.0* 

b  1  2  3.3*  0.36 

3  13  3.35  0.99 

*  3.34  0.96  0.52  0.20 


c  1  2  3.47  0.33 

3  U  3.47  0.34 


3.68  0.52  -0.13  -0.14 
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0.23 
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pattern  apparent  In  the  differences  between  methods;  compare  650  (a) 
where  aethod  1)  and  aethod  4)  agree  well*  to  750  (c)  where  3)  and  4) 
closely  approach  one  another  and  800  (a)  where  1)  and  3)  almost  coin¬ 
cide. 

Thus ,  no  criteria  for  the  beat  fitting  procedure  emerges  from 
the  R  and  a ^  values.  The  additional  amplitudes  a^  and  a^  are  of 
little  assistance,  varying  from  -0.30  to  40.52  at  750  psec.  As  a  re¬ 
sult  the  fitting  procedure  was  chosen  on  a  strictly  practical  basis. 
Method  1)  was  the  simplest  to  employ  and  was  used  in  the  remaining  data 
analysis . 


D.2  Radius  Profile 

In  order  to  generate  R(t)  and  R(t)  values,  as  well  as  Interpolate 
the  discrete  R(t)  data  measured  by  the  procedures  of  the  preceding  sec¬ 
tion  a  smooth  fit  was  made  to  all  the  experimental  data  simultaneously. 

The  shape  suggested  by  the  radius  data  in  Figures  22  and  24-27  is  that  of 
a  sine  wave,  so  a  Fourier  series  was  fit  to  the  makers .  Parameters  which 
were  optimised  through  a  nonlinear  least  squares  analysis  were: 

t  -  the  initial  time  at  which  R(t)  -  0 
o 

tp  -  the  period  of  the  primary  wave 

Cj  -  the  amplitudes  of  the  various  waves 
for  J  ■  1,  2,  . . . ,  5* 

The  functional  form  of  the  wavs  was 


D- 

5  t-t 

Rfit(t)  ■  ^  Cj  sln  ("T£  *  V  0>.3) 

where  are  integers.  The  were  also  varied,  and  the  best  fit 
usually  occurred  for 

*1  "  1 

hm* 

K3  -  3  (D.4) 

K4-5 

S’7 

and  the  contribution  of  these  last  two  waves  was  small.  Here,  best  fit 
means  those  values  for  which  the  sum  of  the  square  residues,  p2,  is 
at  a  minimum 

■>’  -j, >*<'«)  -  W**”’  »'5> 

where  n  is  the  number  of  data  points  and  t£  is  the  time  at  which  these 


data  exist. 


APPENDIX  E:  Additional  Data  for  Photo  Sequences 
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The  camera,  tank  and  laser  are  oriented  so  that  the  laser  enters 
the  fluid  from  the  right  aide  In  each  photograph.  The  bean  is  focussed 
by  the  lens  to  a  "spot",  manifest  as  the  bright  streak  in  each  photo. 
This  streak  is  useful  as  a  reference  point  since  it  appears  at  the 
center  of  the  original  bubble. 

Figure  21. 

Medium  lens  t  1.5  co  from  center 
Laser:  80. OZ,  no  filter 

Figure  22. 

Medium  lens : 

Laser:  80. OZ,  medium  filter 
Figure  23. 

Medium  lens 

Center  to  solid  vail  ^  0.7  1  0.05  cm 
Laser:  97. 0Z,  medium  filter 

Figure  30. 

Distilled  water 

Lens  separation:  8  cm 

Laser:  97. 0Z,  no  filter 

Center  to  solid  vail  0.7  ±  0.05  cm 

Figure  31. 

Distilled  vater 

Lens  separation:  5  cm 

Laser:  97. 0Z,  no  filter 

Center  to  solid  vail  *  0.5  ±  0.1  cm 


